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Abstract 

In this article, we discuss formal invariants of singularly-perturbed linear differential systems in 
neighborhood of turning points and give algorithms which allow their computation. The algo¬ 
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Let d — -j~ x and consider the decisively simple-looking differential equation [34, Introduction, 


Subsection 1.5] 


e 2 d 2 f 2 - (L - e)f = 0 


over the (x, e)-region D : |x| < o-q , 0 < |e| < £ 0 for some positive constants ao and s 0 . Setting 

/ 


F = 


edf 


, this equation can be rewritten as the following first-order system: 


s 8F = A(x,s) F = 


x 2 -e 0 


F. 


(1) 


One can observe that the Iordan form of Aq(x) := A(x, 0) is not stable in any neighborhood of 
x = 0. We refer by turning points to such points where the Iordan form of Ao(x) changes, i.e., 
either the multiplicity of the eigenvalues or the degrees of the elementary divisors change in a 
neighborhood of such points [40, p. 57]. For example, x - 0 might be a turning point for (1). 
Our goal is to compute a formal solution in a neighborhood of x - 0 as e —» 0. 


*A part of this work was developed within the author’s Ph.D. thesis at the University of Limoges, DMI. 
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\x\ = n.Q 



In the region /?o|e|^ 3 < |x| < ao where /3o is a positive constant, the transformation 

transforms system (1) into 


F -T\G where T\ = 


0 x 3 ^ 2 


(x^e)x 3/2 dG = { 


0 1 
1 0 


+ (x 3 e) 


-1 4x 1/2 


)G, 


which can be rewritten, by a self-explanatory change of notation x 3 e = 4 ; , as 


gx 3/2 dG = { 


0 

1 


0 0 

1 

0 

-1 fx 1 ' 2 


A transformation G — TnU where Ti = 
diagonalized system 

tx y2 dU = { 


1 i 

2 2 
i T 

'2 2 J 


-1 0 

0 1 


}G. 

+ 0(1;) would then result in the block- 

+ omu. 


Consequently, the system splits into two first-order linear differential scalar equations. We 
can then write down the following fundamental matrix of formal solutions for (1) in this 
region: 


Fouter = TiT 2 exp( 


-fx 5 % 1 +0(x- 1 ' 2 ) 

0 


2 r 5/2„-l 


0 

+ 0(x~ 1/2 ) 


• In the region |x| < Pq\s\ 1 ^, we shall perform a so-called stretching transformation. Namely 
we set r = xs ~ 1/3 and d T = Then, for all r such as |r| < oo, except possibly for 
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neighborhoods of the roots of r 2 3 - 1 = 0 2 , the transformation F - L\ G where 


U = 


reduces system (1) to 


1 0 

0 e l ' 2 


0 


1 


r 3 - 1 0 


G. 


s ll6 d T G = 

And the transformation G - LoU where 

\ + i r 3 + 0(t 6 ) \ | r 3 + 0(t 6 ) 

-{ + ^ 3 + 0(t 6 ) 5 - |r 3 + 0(t 6 ) 


L2 = 


reduces system (2) to 


1/6 d T U = { 


-i + ^r 3 + 0(t 6 ) 


0 


+ 0(e l/6 )}U. 


0 i - kj s + 0(T b ) 

A fundamental matrix of formal solution of system (1) is then given by: 


= L\ L 2 exp( 


-;r+(i/8)T 4 +0(T 7 ) n 

K l/f, u 

iT-(il&)T A +Q(T 1 ) 


0 


.1/6 


)■ 


( 2 ) 


We call F ouler (resp. F, rmer ) an outer (resp. inner) solution as it is obtained in the outer (resp. 
inner) region around x = 0. The corresponding differential systems are sometimes referred to 
as outer and inner differential systems as well. However, unlike the above particular example, 
intermediate regions might be encountered. 

In this article, we study the construction of formal solutions (the formal reduction), of sys¬ 
tems of the following form: 

e h dF = A(x, s)F, (3) 

where the entries of the matrix A(x, e) are formal power series in e, whose coefficients are formal 
power series in x\ and h is an integer which we call the e-rank of the system. 

In [? , p. 74], Iwano summarized as follows the problems needed to be resolved to obtain 
the complete knowledge about the asymptotic behavior of the solutions of a singularly-perturbed 
linear differential system: 

(PI) Divide a domain [D] in (x, e)-space into a finite number of sub-domains so that the solution 
behaves quite differently as e tends to zero in each of these sub-domains; 

(P2) Find out a complete set of asymptotic expressions of independent solutions in each of these 
sub-domains', 

(P3) Determine the so-called connection formula; i.e. a relation connecting two different com¬ 
plete sets of the asymptotic expressions obtained in (2). 


2 The turning points other than r = 0, namely the roots of this equation do not explicitly correspond to the original 

turning point. They are referred to as secondary turning points [34], 
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In order to resolve (PI) and (P2), Iwano proceeded by associating a convex polygon to first- 
order systems in analogy with the scalar case, after imposing a precise triangular structure on 
A(x, s). (P3) remains generally unresolved and it is out of the scope of this article as well (see, 
e.g. [20, 22, 38] and references therein). The hypotheses on A(x, e ) were eventually relaxed 
for special types of systems in a series of papers by Iwano, Sibuya, Wasow, Nakano, Nishimoto, 
and others (see references in [40]). A famous and prevalent method among scientists is the Wkb 
method (see, e.g. [26]), which, roughly speaking, reduces the system at hand into one whose 
asymptotic behavior is known in the literature. However, as pointed out in [39], in view of the 
great variety and complexity of systems given by (3), it is rarely expected that a system can be 
reduced into an already investigated simpler form. 

The methods proposed in the literature for the symbolic resolution of such systems are not 
yet fully algorithmic, exclude turning points (see [18] and references therein), treat systems of 
dimension two only, rely on Arnold-Wasow form [3], and/or impose further restrictions on the 
structure of the input matrix. Moreover, there exists no computer algebra package dedicated 
to the symbolic resolution of neither system (3) nor the scalar n^'-order scalar equation. The 
widely-used softwares Maple and Mathematica content themselves to the computation of outer 
formal solutions for scalar equations, and so does [27], 

In this article, we attempt to resolve (P2) algorithmically, that is give an algorithm to construct 
a fundamental matrix of formal solutions of an input system. We also attempt to give an insight 
into (PI) within some examples. 

As illustrated by the introductory example, the formal reduction of a singularly-perturbed 
system given by (3) leads inevitably to the consideration of more general systems. These consid¬ 
erations will be discussed in Section 1. In Section 2, we give necessary preliminaries including 
results on n ,h -order scalar equations, the base cases, and the well-known Splitting lemma which 
splits the system into subsystems of lower dimensions, whenever possible. We discuss our con¬ 
tribution in the remaining sections, as we deviate from the classical reduction to give algorithms 
which act on the system directly without resorting to an equivalent scalar equation or to the 
Arnold form: First we refine in Sections 3.1 and 3.2 the algorithms which we proposed in [2] for 
the resolution of turning points and reduction of the s-rank h to its minimal integer value. Then 
in Section 4, we associate an e-Newton polygon to a given system and show its role in retrieving 
formal invariants. And finally, we make use in Section 5 of these notions introduced to propose 
an overall formal reduction algorithm. 

As we have already indicated, the literature on this problem is vast. In the hope of keeping 
this article concise but self-contained, we restrict our presentation to constructive results which 
contribute directly to our process of formal reduction. The omitted proofs are collected in the 
appendix along with illustrating examples. 

In this paper, we content ourselves to the process of formal reduction. Any reference to the 
asymptotic interpretation of formal solutions will be dropped in the sequel. One may consult in 
this direction the method of matched asymptotic expansions (see, e.g. [37, 38] and references 
therein) or composite asymptotic expansions (see, e.g. [20, 22] and references therein). 

In the sequel, we adopt C as the base field for the simplicity of the presentation. However, any 
other commutative field T of characteristic zero (Qcf c Q can be considered instead. The 
algorithms presented herein require naturally algebraic extensions of the base field. However, 
such extensions can be restricted as described in [8, Section 7.2], 

Since we are leading a local investigation, we assume that the input system has at most 
one singular point, otherwise the region of study can be shrunk. Moreover, we place this singular 
point at the origin, otherwise a translation in the independent variable can be performed (x —> 1 /x 
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for a singularity at oo). 


Notations 

• C is the field of complex numbers, Q is the field of rational numbers, and Q _ - {q e 

Q \ q < 0). 

• We denote the dimension of the systems under consideration by n. We use r for the alge¬ 
braic rank of specified matrices (leading matrix coefficients) and h for the e-rank. 

• We use upper case letters for algebraic structures, matrices, and vectors; the lower case 
letters i,j,k,l for indices; and some lower case letters (greek and latin) locally in proofs 
and sections, such as w, v, cr, p, g, etc... 

• We use x, t, r as independent variables and d to signify a derivation. Since x dominates this 
article, we drop it in the derivation. We set d — and d T = 4^ . 

• We use F for the unknown «-dimensional column vector (or n X n-matrix); we use G, U, 
W, and Z in this context as well. For scalar equations, we use / and g. 

• Idxn (resp. /„) stands for the identity matrix of dimensions d x n (resp. n X n); and 0,j x „ 
(resp. ()„) stands for the zero matrix of dimensions d x n (resp. n X n). If confusion is 
unlikely to arise, we simply use 0 to denote Odxn or 

• We say that M e A1„(R) whenever the matrix M is a square matrix of size n whose entries 
lie in a ring R. 

• GL„(R) is the general linear group of degree n over R (the set of n X n invertible matrices 
together with the operation of matrix multiplication). 

• C[[x]] is the ring of formal power series in x whose coefficients lie in C; and C((x)) is its 
fraction field. The usual valuation of an element g(x) of C((x)), i.e. the order of g in x, is 
denoted by val x (g), with val x { 0) = +°o. The ring C[[x]] (resp. the field C((x))) is endowed 
with a differential ring (resp. field) structure by considering the derivation d. 

• C[[e]] is the ring of formal power series in s whose coefficients lie in C; and C((e)) is 
its fraction field. C[[x]][[s]] (resp. C[[x]]((e))) is the ring of formal (resp. meromorphic) 
power series in s whose coefficients lie in C[[x]]. Similarly, C((x))[[s]] (resp. C((x))((e))) 
is the ring (resp. field) of formal (resp. meromorphic) power series in s whose coefficients 
lie in C((x)). 

• Any element / e C((x))((e)) can be written as / = tezA(x)e*, where fk(x) G C((x)) 

for all 1 6 Z. The derivation d extends naturally to C((x))((e)) by the formula df — 
2fez(d/)t(x))s A '. Moreover, the map val e : C((x))((e)) —> Z U oo defines a valuation over 
C((x))((s)), satisfying the following properties for all /(x, e), g(x, s) in C((x))((e)): 

- val E (/) = oo if and only if / = 0; 

- val e ( fg ) = vale (/) + vale C?); 

- vale (f + g) > min ( val e (/), val E (g)). 
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• We give the blocks of a matrix M with upper indices, e.g. 

|M U M 12 

M - [ m21 m 22 J, 

and the size of the different blocks is dropped unless it is not clear from the context. 

1. Ring of coefficients 

In this first part, we investigate a suitable ring of coefficients to treat first-order singularly- 
perturbed linear differential systems and n th -order scalar singularly-perturbed differential equa¬ 
tions, in a neighborhood of a turning point (see also [40, Chapter 2] and references therein). For 
( Uq, Vq) e R 2 , we put 


P(U 0 , V Q ) = {(£/, V) e R 2 | U > U 0 and V > V 0 ). 

Let / = Etez 6 C((x))((e)) and let Pf be the union of P{k,val x (fi)),k e Z. Then the 

Newton polygon of/, denoted by Nf, is the boundary of the convex hull in R 2 of the set Pf. We 
now consider, for cr, p e Q with cr < 0, the half-plane 

H^ p = {(U, V) e R 2 | V > crU + p}. 


One can verify that 

7C = { f e C((x))((e)) | Nf c H (r , p for some cr e Q",p e Q } 

is a differential field endowed with the derivation d — d/dx and it is the field of fractions of the 
ring V. - 7C n C((x))[[e]]. Moreover, the elements of 7C can be characterized geometrically as 
follows: 

/ = Yj e7C^3creQ-,peQ|Vk, val x (f k ) > crk + p. 

k£ Z 

In fact, given a non-zero element / of 7C, we set v/ := val E (f). Among the half-planes which 
contain Nf, we consider the half-planes which are bounded below by a straight line passing 
through the point (v/, val x (f V[ )), i.e. the half-planes H ljpf with pf := val x (f Vf ) - crvf. Then, 
among the H lrp \, we pick the half-plane of maximal slope, which we denote by cry. Evidently, 
H(T f . P! is determined by a straight line which passes through the point (vj, val x (f Vf )) and point(s) 
(k , valjfk)) for at least one k > v. We denote by (Sf) this straight line whose equation is given 
by V = (TfU + pf, and we say that it is associated to /. For f - 0, we set cry = 0 and pf = 0. 


6 





* a 

x g 

3 f 

0 s 

- V—3U 

-V—11U+14 

v -5i; : 

- V—3U-2 



Figure 2: Geometric interpretation of Example 1: (S a ), (So). (Sf), (S s ). are traced in black, green, blue, and red 
respectively. 

Example 1. Given a nonzero element f of C((x))((e)). If we can construct a straight line of 
finite slope which passes through (Vf, val x (f Vf )), and stays below the points (k,val x (fk)) for all 
k > v, then f e'K. We illustrate the following examples in Figure 1: 

• Let a = x 3 e + x~ 8 s 2 = x 3 e(l + x _11 e) then cr a = -11. 

• Let g = YikLo then cr g = -3. 

• Let f — x 2 + x~ 3k ^ then (Tf — -5. 

• Let s = x~ 2 + x~ 3k s^ then cr s = -3. 

The following two lemmas give further insight into 7<7 

Lemma 1. If f(x, £) e ‘K is non-zero, then there exists e e Z~ such that /(x, x e s ) e C[[x]]((e)). 

Proof. Let / = Yff =v , fk(x)£ k . We set pk = val fff) and fi = xP k fk. Then for any e e with 
e < cr we have: 

/(x,e) = s Vf (f Vf + f Vf+1 S +f Vf+2 £ 2 + fy f+ 3 S 3 + ...) 

= £ Vf xf Jr f <J Vf + x Mv / +1_/ ' v // V/+1 e + x Mv f +2 ~ ll '’f f Vf+ 2 S 2 + x^ v / +3 " / ' v // V/+3 e 3 + ...) 

= £ Vf X Mrf (f Yf + xf' , f* 1 ~t‘ y r e f Vf+ fx e £) + xf ,v f* Z ~t‘' , r 2e f Vf+ 2 (x e S ) 2 + tf v P 3 ~ tIr f~ 3 e f Vf+ 3(x e £) 3 + ...). 

But, pk > CTf{k - Vf) + p Vf for all k > Vf. Hence, pk > e(k - vf) + p Vf for all k > v/. We can thus 
choose e to be the largest integer which is less than or equal to <jf. □ 

Lemma 2. The field of constants of the differential field PK, d) is C((e)). 

Proof Let / e 7f and suppose that df = 0. Then, Yff =v (dffix'fif = 0 which yields dffix) = 

0 for all k > Vf. Thus, for every k > v/, fifix) = q for some q e C, which yields / = 

■ □ 
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Henceforth, algorithmically, we can use the following handy representation of any element 

/ c '/< : 

oo oo oo 

fix, e) = f k (x)s k = ^ x~ k(r 'f k (x)(x (ri £) k = x Pl ^ x~ kai ~ Pl f k (x)f k f , with f f = x^s. 

k—Vf k—Vf 

(4) 


=V/ 


Evidently, val x (x (r,k Pf f k (x)) > 0 for all k > v/. Moreover, under the notation ff, we have: 


- Yj (d( f k) + -j-fiXf- 

k—Vf k= v f 


1.1. Input systems and input equations 

Motivated by the introductory example, and with the help of the above notations, we will 
investigate the following system rather than (3): 


dF = A(x, e)F where A e Mni'K). 


(5) 


The valuation val e (and the usual valuation val x of C((x))) extends naturally to matrices A — 
(ai j ) e Mni'K) by val e (A) = min i ; val e (aij). We can then define the Newton polygon of A e 
Mni'K), va e Z, pa 6 Q, and cta <£ Q~, in an analogous manner to those of / e 'K. Hence, we 
can express A(x, s) as follows: 


A = X AkE A > A *a * °, Pa ~ valxiAyf) - o- A v A , and val x (A k ) > cr A k + p A V k> v A . (6) 

k=v A 

Hence by setting = xF A s, we can write in analogy with (4): 

OO 

A(x,s) = x PA X x- krrA - pA A k (x) 

k=v A 

OO 

= Y i x ~ krA ~ vaUA,A) Av A 4*{x)& where x~ valM ^A Va \ x=0 * 0. 

k =0 

Thus, we can rewrite system (5) as follows: 

oo 

C A x~ PA dF = (J] x ~ krTA _vai ' (/t '’A >a va +k (x) f k A ) F, where x~^- valAA ^A VA+k e M n {C[[x}}). 
k =0 

Similarly, if we consider a n rh -order differential equation whose coefficients are elements of 7C: 

a n (x,s)d n f + a„-i(x, e) d" -1 /+... + a | (x, e) df + a 0 (x,s) f = 0, (7) 

then setting F - [/, 3/, <9 2 /,... ,d n ~ 1 f] T , we can express (7) equivalently as a system <7F = 
A(x,e)F where Aix, s) e Mni'K ) is a companion matrix, and we define v a v A , p a := /r,i, 

<r a ■= Oa, and := % A . 
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1.1.1. Notations 

Thus, in this paper, we treat system (5) in the following form: 

oo 

[Act J ? a x Pa 8F = A(x, £ A )F A(x, £ a ) = 2 A k {x)£ k A g M„(K), A 0 (x = 0) * 0, 

fc =0 

( 8 ) 

where ^ = x <TA e , cr^ g Q~, p A G Q, and h G Z. Under the definition of f A , A(x,£ a ) = 
Z k=0 A k (x)& with A k (x) G M n ( C[[x]]) for all it > 0. We refer to Ao(x) as the leading coefficient 
matrix and to Ao.o := A a (x = 0) as the leading constant matrix. 

We also consider n th -order scalar linear differential equations of the form: 

[« C J &'f + a n -i(*,£,)d" -1 / + ... + ai(x,^ a )df + a 0 (x,£ a )f = 0, (9) 

where fl,(x, £ a ) = a i,k(x) e* G 7f for all i G {0,...,«}, with G C[[x]] for all k >0. 

We remark that a system [/f fr J given by (8) can be also expressed as a scalar equation of 
the form (9). The theoretical possibility stems from the work of [19] which discusses cyclic 
vectors. Moreover, an algorithm to compute a companion block diagonal form of [A iri ] (and 
hence equivalent scalar form) can be obtained by generalizing the work in [7] developed for 
unperturbed systems (see Appendix A). 

In the sequel, for the clarity of the presentation, the index A (resp. a) will be dropped from 
v A , Pa, o~a, and g A (resp. v a , p„, cr a , and f„) whenever confusion is unlikely to arise. 


1.2. Equivalent systems 

Consider system [A^J given by ( 8 ). Let T g GL„( r K ) then the transformation F — TG (also 
called gauge transformation) yields a system 

[A cr .] Zz h x p dG = A(x,f A )G, where A e MM). (10) 

for some cr A G Q“, p G Q, and h G Z. We say that systems [A^J and ^[A^] := [A^] 
are equivalent. Examples of such transformations and their applications are the transformations 
T \, To, L\ , Lo applied within the introductory example. 

In the sequel, we seek at each step a transformation which yields for an input system [A rTA ], 
an equivalent system [A^], so that either h < h or [A, r .] can be decoupled into systems of lower 
dimensions. Using a recursive process which employs this approach, we aim to break down any 
input system [A^J into system(s) which can be treated with known methods (see basic cases of 
Section 2.1). 

Evidently, cry. might differ from <:r A and hence £4 will be updated to £y after the application 
of T (see the introductory example), which brings us to the next subsection. 

1.3. Restraining index and (PI) 

As we are lead naturally to the study of system [A^-J with cr A g Q~ within the study of 
system ( 8 ), we will need to understand the inevitable growth of the order of the poles in x within 
the reduction. For instance, in the introductory example, we start with an input system for which 
cr a = 0 and arrive at an decoupled system with cru — -3. 

As described in [40, Chapter 2] and as we have observed, these poles grow at worst linearly 
with k. From here initially stems our main motivation for our choice of the ring of coefficients 
and the representation in terms of £ A - the orders of the poles in x which might be introduced in 

9 


an input system within the process of reduction, are stored in p and cr, which allows investigating 
them at any step within the reduction process. We can thus talk about a restraining index p = 
-1 /cr. This is the first step towards determining the inner and intermediate regions. In fact, the 
complicated behavior anticipated in the neighborhood of turning points can be investigated with 
the help of a sequence of positive rational numbers: 

[p] 0 = po < Pi < P 2 < ■ ■ ■ < Pm- (11) 

With this sequence, the domain \x\ < xo of [D] can be divided into a finite number of sub-domains 
in each of which the solution of behaves quite differently (see [23, Intro, pp 2] and [? , Chapter 
1]). In case [p] is known, one can apply adequate stretching transformations to the original 
input system, i.e. a change of the independent variable of the form r = xs~ Pi , i e [ 1 ,..., /??). 
Then, using our proposed algorithm, we can construct inner and intermediate solutions. In this 
paper, we do not investigate (PI), i.e. we do not compute [p]. However, we motivate a plausible 
approach in some examples (see Appendix E). 

2. Preliminaries 

A detailed discussion of the scalar case (n = 1) and regularly-perturbed systems (/? < 0), is 
given in [18, Chapter 3, Section 1, p. 52 -56] for cr A = p A = 0. In this section, we give a brief 
discussion which discards this restriction on cr A and p A . 

2.1. The base cases 

2.1.1. The case h < 0 

If /; < 0 then the solution of system [/\, r J can be sought, up to any order p, upon presenting 
the solution as a power series in = xF A s, i.e. 

CO 

F = Yj F k{x) %A h+k - 
k =0 

The latter can then be inserted into [A (rA ] £ l x p dF = AF and the like powers of £4 equated. This 
reduces the symbolic resolution to solving successively a set of inhomogeneous linear singular 
differential systems in x solely: 

x p 8F 0 = (A 0 (x) + h cr a x p ~ l I) F 0 

x p dF\ = (Ao(x) + (h - 1) cr a x p ~ l I) F\ + A\(x) Fq 

x p dF v = (A 0 (x) + (h - p) cr a x p ~ l I) F M + Ai(x) + ... + A M (x) F 0 . 


For any k > 0, the solution of the homogeneous system singular in x can be sought via the 
Maple package Isolde [14] or the Mathemagix package Lindalg [29] (which are both based 
on [8]). Afterwards, the solution of the inhomogeneous system can be obtained by the method 
of variation of constants (see, e.g. [5, Theorem 3, p. 11]). We remark that a transformation 
Fk — a ka ' G reduces the inhomogeneous system in the dependent vector /•’/. to: 


x p dG = (Ao(x) + h a -a x p 1 /„) G(x ) , 
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which coincides with that of Fq(x). Thus we have: 


F k (x) = x~ k<TA Fo(x) J Fg 1 (x)x kcrA [Ai(x) F k .fx) + ... + A k (x) F 0 (x)] dx. 

2.1.2. The scalar case 

The n th -order singularly-perturbed scalar differential equations are treated in [27] upon in¬ 
troducing an analog to the Newton polygon and polynomials (see e.g. [6] and references therein 
for the unperturbed counterparts of these equations): the e-polygon and e-polynomials. In this 
subsection, we adapt this treatment to the more general equation [a cr fi given by (9): 

[flo-J d"f + a„-\(x,f a )d n ^f + ... + a\(x,^ a )df + a 0 (x,fi a )f = 0. 

We first give an analog to [27, Definition 2.1] and then generalize [27, Proposition 4.1], 
whose proof is easily adaptable from [27] (see Appendix B). 

Definition 1. Consider the scalar equation [a o-J given by (9). For (uq, Vo) e R 2 , we put 
P(U 0 , Vo) = {((/, V ) € R 2 | U < U 0 and V > V 0 ). 

Set V,- = val E (ai) for i e [0,..., n) and let P a be the union of P(i, v,-). Then the e-polygon of [flo-J, 
denoted by N e (a), is the intersection of R + X R with the convex hull in R 2 of the set P a . We denote 
the slopes of the edges ofN e (a ) by [ei,..., ei). These slopes are non-negative rational numbers. 
And, for every j € [1 ,...,(}, consider the algebraic equation given by 

e 

t Ej ) J^a^^^ix) = 0, 

k=0 

where 0 < ;’o < i\ < ■ • ■ < it — n denote the integers ifor which ( i , v,) lie on the edge of slope ej 
of the e-polygon, and a iv .(x) = f~ Vi afx, f a )\e a = 0 - W ’e say that (Ej) is the e-polynomial associated 
to the slope ej. 

Proposition 1. Consider a nonzero [a, r f given by (9) and its e-polygon N r (a) of slopes {e \,..., ef. 
Let 

f(x,£ a ) = exp( f q(x^a)dx) with q(x^ a ) ± 0 € [J C((x))((^ /5 ))- 

^ seW 

If f(x,f a ) satisfies [a^j] formally then 

q(x,fa) = -JjWx) + O(ja)) 

ha 

for some j € [1 ,...,(}, and X(x) is one of the non zero roots of the e-polynomial (Ej) associated 
to ej . 

The full expansion of q(x,f a ) can be obtained with successive substitutions of the form / = 
exp(f -dx) g [27, Proposition 4.2]. 
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Figure 3: The e-polygon associated to (12) in Example 2 


Example 2. Consider the following linear differential equation with cr a = 0: 

M d 2 f - fd 2 f - = o. (12) 

We wish to study this equation according to its e-polygon (see Figure 3). We have two slopes: 

• The slope e\ = | for which (E\) x X 2 + 1 = 0 whose nonzero solutions are X = ±-4=. 

• The slope e 2 = 2 for which (Ef) X 3 - xX 2 = 0 whose nonzero solution is X — x. 

Thus, the leading terms (of the exponential part) of solutions of ( 12) are given by exp (f pdx) = 
exp (f dx) and exp (+ f }/1 dx) = exp (+ J ^ 3/2 dx). 

One can also refer to [27, Section 5] for examples on the harmonic oscillator with damping 
and Orr-Sommerfeld equation. 

Based on the above, we give the following definition: 

Definition 2. Consider the scalar equation [t; fr J given by (9). We call the largest slope of the 
e-polygon N E (a), the e-formal exponential order (e-exponential order, in short): 


u> £ (a) = max {ej, 1 < j < l}. 


Clearly, under the notations of this subsection, we have: 


n— 1 

ajJa ) = max ( 0 , - 

i=0 


val e (ai) 
n - i 


2.2. Consequences of the equivalence between a scalar equation and a system 

Due to their equivalence, the formal solutions of a first-order system [A^] given by ( 8 ) can 
be computed from an equivalent scalar equation [a, r J given by (9), using the e-polygon of [a, r J. 
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Moreover, (11) can be computed as well using the Iwano-Sibuya’s polygon (which collects valu¬ 
ations with respect to both x and e [23]). However such a treatment is unsatisfactory (see e.g. [27, 
Conclusion]) since it overlooks the information that we can derive from the system directly and 
demands an indirect treatment. Nevertheless, it plays a key role in the theoretical basis of the 
algorithm which we will develop in later parts. Most importantly, the equivalence between [,4, r J 
and [a,,.], and the invariance of the e-exponential order under gauge transformations, allows the 
definition of the e-exponential order of system [A fr J as follows: 

Definition 3. Consider a system [A^] given by (8) and an equivalent scalar equation [a^-J given 
by (9). The e-polynomials and s-formal exponential order a) e (A) o/tA^] are those of[ao-J. 

It follows that, given a system [A aA ], we have: 

—vaL (A,-) 

w fi (A) = max (0, -(13) 

0<i<n n — l 

where [a^J is a scalar equation equivalent to [A^J. 

In this paper, we give a direct treatment of the system which will lead eventually to the com¬ 
putation of (j) t: {A ), the e-polynomials, and consequently a basis of the space of formal solutions, 
without resorting to the equivalence of an input system to a scalar equation. Our treatment relies 
on the properties of the eigenvalues of the leading coefficient of the input system. We distinguish 
between three cases (distinct, unique, or zero eigenvalues), the first of which is classical and is 
recalled in the next subsection. 


2.3. e-Block Diagonalization 

Consider a system [A, Fa \ given by (8). A classical tool in the perturbation theory is the 
so-called splitting which separates off the existing distinct coalescence patterns. Whenever the 
leading constant matrix Ao,o := Aq(x = 0) admits at least two distinct eigenvalues, the system 
can be decoupled into subsystems of lower dimensions: 

Theorem 1. Consider system [A fr J given by (8) 

[A,rJ t h A x P dF = A(x,f A )F, A(x,^)gM„(7?), A 0 ,o * 0. 


If A 0 ,o 


4 11 
^ 0,0 

O 


O 

A 22 

A 0,0 


such that Apl and A^ have no eigenvalues in common, then there exists a 


* 0,0 


unique transformation T(x,£t) = E^=o Fk(x)fj € GL n ( r R), given by 


k=0 


r(*,fr) = / + 2 

such that the transformation F — TG gives 
[Ao-j] % l ~x p dG = A(x,Ta)G = 


O T l k \x) 
Tf{x) O 


fy k . 


A n (x,&) 

O 


O 

A 22 (x,&) 


G where Aq.o = Aq,o- 


Moreover cr T , cr^> cr where cr ■ 


( <r a + V (/><! 
I cta otherwise 
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Well-known proofs proceed in two steps: Block-diagonalizing A 0 (x), followed by successive 
block-diagonalization of coefficients (see, e.g., [21, Theorem XII-4-1, p. 381] and [18, Chapter 
3, p. 56-59] for the particular case of cr A — p A = 0, and [40, Theorem 8.1, p. 70-81] for the 
general case). We give a constructive proof in Appendix B. We remark that the specified form 
of Ad i, in Theorem 1 is non-restrictive. It suffices that Ao.o has at least two distinct eigenvalues. 
Such a block-diagonal form can be then attained by Jordan constant transformation. In case the 
base field is not algebraically closed, the weaker transformation given in [15, Lemma A.l] can 
be applied instead, since it does not require any algebraic field extensions. 

After splitting the input system, we proceed in our reduction for each of the decoupled sub¬ 
systems in parallel. We can thus assume, without any loss of generality, that the leading constant 
matrix of the input system [A (Ta | has a unique eigenvalue y 6 C. Thus, upon applying the so- 
called eigenvalue shifting, i.e. 



(14) 


one can verify that the resulting system L[A, r J has a nilpotent leading constant matrix. Thus, it 
remains to discuss this case of nilpotency. This is the goal of the two following sections. 

3. Two-fold rank reduction 

Without loss of generality, we can now assume that system [A^] given by ( 8 ) is such that 
Ao,o is nilpotent. At this stage of reduction, our approach diverges from the classical indirect 
ones which require Arnold-Wasow form, a cyclic vector, or reduction to a companion form. Two 
cases arise: 

• Ao.o is nilpotent but Aq(x) is not: In this case, the eigenvalues of Aq(x) might coalesce in 
some neighborhood of x = 0, which would cause the so-called turning point. In Subsec¬ 
tion 3.1, we propose an algorithm to compute a transformation which reduces the input 
system to a system for which both Ao.o and Ao(x) are nilpotent. 

• Both Ao.o and Ao(A) are nilpotent: In this case, we propose in Subsection 3.2 an algorithm 
to compute a transformation which can reduce the e-rank h of system [A, T J to its minimal 
nonzero integer value. This minimality is an intrinsic property of the system and its so¬ 
lutions. Its computation paves the way for the retrieval of the s-exponential order of the 
system, as we show later in Section 4. 

3.1. Resolution of turning points 

We illustrate the source of turning points in the following example: 

Example 3. [38, p. 223] Let a(x), b(x) be holomorphic for x in a region LI e C and consider 


M(x) = 0 
0 


x 


1 a(x) 
x b(x) 
0 0 
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a(x) f3(x) y(x)(—a(x)x~ l + b( x)x -2 ) 

0 a(x ) -y(x)b{x)x~' 

0 0 y(x) 

functions a(x),[3(x), y(x) such that a(x) + 0 and y(x) + 0 for all x e Cl, yields for x 4 0: 


The transition matrix T(x) = 


, with arbitrary scalar 


J(x) = T~ l M{x)T{x ) = 


x 1 0 
0x0 
0 0 0 


If a(x),/3(x),y(x) are holomorphic, so is T, provided that zero is not a point of Cl. If 0 € Cl, then 
T(x ) has a pole at x = 0 unless the following conditions are satisfied: 

b(0 ) = 0 and db(0 ) = a(0 ). 

Therefore, M(x ) is in general not holomorphic ally similar to 7(x) in regions that contain x = 0. 
Moreover, if the second condition is not satisfied, two possibilities arise: 


ifb( 0) = 0 then there exists a constant invertible matrix S such that S 1 M(0)S - 


0 1 0 
0 0 0 
0 0 0 


Then M(x) is pointwise similar to J{x), even in regions containing x = 0. 

so that J(x) is not holomor- 


ifb( 0) 4 0 then the Jordan form of M(0) is 7(0) = 
phic at x = 0. 


0 1 0 
0 0 1 
0 0 0 


Hence, it might occur that the Jordan matrix is not holomorphic in some region or the Jordan 
matrix is itself holomorphic but, nevertheless, not holomorphically (although perhaps point-wise) 
similar to M(x). Below is an example which exhibits a system whose leading coefficient matrix 
has such complications: 

Example 4. [40, p. 57] Consider the following system with cr A — 0, i.e. £4 = s: 

[ArJ ? A 9F = 


0 

1 

0 


0 

1 

0 

0 

0 

1 

F where Aq(x) = 

0 

0 

1 

U 

0 

X 


0 

0 

X 


It follows from the form of Aq(x) that the origin is possibly a turning point. We remark that this 
system is equivalent to (12), as it is obtained from the latter by setting F = [f,f 2 A df,fi i A d 2 f] T . 

We recall the following proposition which we first gave in [2], and we refine its proof: 

Proposition 2. [2, Section 3] Consider system (8) given by 

00 

[AcrJ f h A x‘>dF = Mx,f A ) F = Y,A k {x)£ k A F 

k=0 

with f A — x ITa e. Suppose that the leading constant matrix A 0 0 is nilpotent and Aq(x ) has at least 
one nonzero eigenvalue. Then there exists an invertible polynomial transformation T in a root of 
x such that the transformation F — TG and a re-adjustment of the independent variable x results 
in a system whose leading constant matrix has at least one nonzero constant eigenvalue. 
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Example 5. Consider the matricial form of Weber’s equation (o~a = 0, p A = 0) given by: 




y dF_ _ 
^ dx 


0 1 

x 2 0 


F, 


with = xF A e. Let T(x) - 


1 0 
0 x 


then F — TG yields: 


dG 

[B - ] % = 


X -ftx' 


.-1 


G. 


Upon factorizing x, we have: 


e -\ dG 

ZaX — = 
ax 


0 1 

1 -ftx -2 


G. 


We compute crj; = -2 and thus set £4 = x ~e. The above system can be rewritten as: 




hA X 1 ~ 

dx 


0 1 
1 -ft 


G. 


Clearly, Ao,o has 2 constant eigenvalues and so the splitting lemma can be applied to decouple 
this system into two scalar equations. 

Proof (Proposition 2) The eigenvalues of Ao(x) admit a formal expansion in the fractional pow¬ 
ers of x in the neighborhood of x = 0 (see, e.g., [25]). We are interested only in their first nonzero 
terms. Let p(x) = JfJ=oPj x ^ s be a nonzero eigenvalue of Ao(x) with 5 e N*, jAs = 1, and whose 
leading exponent, i.e. smallest j/s for which pj + 0 , is minimal among the other nonzero eigen¬ 
values. Without loss of generality, we can assume that s = 1, otherwise we set x — t s . By [9, 24], 
there exists T e GL„(C(x)) such that for Bq(x ) = T~ 1 Aq(x)T = Bo.o + Bi,o* + ...: 

0b o(x) (A) = x' ank(Boo) det(/17 + ^ + Bi,o)L=o 

x 

does not vanish identically in A. Let v > 0 denote the valuation of Bq(x) in x. By [24, Proposition 
1], there are n - deg(P) eigenvalues of B<fx) whose leading exponents lie in [v, v + 1 [, and deg(0) 
eigenvalues for which the leading exponent is equal to or greater than v + 1. Then, applying the 
transformation F = T(x)G to the system [4, r J yields 


[«,rj t A x P dG = B(x,{ a )G = Yj B k(x^ k A G 

k= 0 

where Bq(x) := x v (B qq + B 1 qx +...) and B 0,0 has n-deg(P) eigenvalues whose leading exponents 
lie in [0,1[. Let A(x, = x _y B(x,£ 4 ) then the order of poles introduced in x in A(x, £ 4 ) is at 
worst (span v (r)+v), where span t (T) is the difference between the valuation and the degree of the 
polynomial transformation T(x) in x. Hence, cr^ > <r A - span t (T) - v and | B ITa ] can be rewritten 
as 

OO 

[Ar A ] fy? 8G = A(x, ft) G = YjA k (x)£ k A G 

k= 0 
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where £4 = x°~ A s, peZ, A^(x) e A1„([[x]]), Ao(x) = x~ v Bu(x), and Ao,o = # 0 , 0 • Then p(x) is an 
eigenvalue of A {) (x) with a minimal leading exponent and hence it is among those whose leading 
exponents lie in [0,1[. By our assumption s = 1, and hence the leading exponent of p(x) is zero 
and /io 4 0. Since jiq is a nonzero eigenvalue of Aq.o, it follows that the Ao.o is non-nilpotent. □ 

We remark that the eigenvalues of Ao(x) are the roots of the algebraic scalar equation fix. A) = 
det(Ao(x) - AI„) = 0 and can be computed by Newton-Puiseux algorithm. The sought polynomial 
transformation T can be computed via Isolde, miniIsolde, or Lindalg. 


Remark 1. Proposition 2 leads to the following observation about the detection of a turning 
point: Consider an input system [A (TA ] given by (3) with <t a — 0, Pa — 0, and restraining index 
Pa = -1 /o' a- If [Act a ] has a turning point at x = 0 then, by the end of the process of formal 
reduction of this system (i.e. whenever all decoupled subsystems are either of non-positive s- 
rank or are scalar equations), at least one of its decoupled subsystems has a finite restraining 
index. 


Example 6 . Consider the following system whose (Ta - 0 and if — £■' 


[^J 


QdF = 


0 

1 

0 


0 

1 

0 

0 

0 

1 

F where Aq(x) = 

0 

0 

1 

iA 

X 

0 


0 

X 

0 


We first compute s = 2. We set x = t z and compute T = 


T = 


0 0 1 

0 x 1/2 0 

0 0 x 


. Then F = TG yields: 


0 0 
0 t 
0 0 


Or equivalently, we consider 


i A dG = x 


- A/ 2 j 


0 

1 

0 


0 

0 

0 


0 

0 

0 

1 

0 

0 

+ 

0 

0 

0 


0 

-1 

0 

0 

1 

0 


1 

0 

0 


0 

0 

-2 


G , 


which we rewrite as: 





0 

1 

0 


0 

0 

0 


0 

0 

0 

dG ={ 

1 

0 

0 

+ 

0 

0 

0 

iA + 

0 

-x 2 

0 


0 

1 

0 


1 

0 

0 


0 

0 

- 2 x 2 


d ) G, 


where £4 = x 3 ^ 2 e. Now that the leading term is a constant matrix with three distinct eigenvalues, 
we can proceed by applying the Splitting of Theorem 1. 

In the above example, the system could be decoupled thoroughly into three scalar equations. 
However, we might encounter different scenarios as well, as illustrated in the following example. 


Example 7. [40, p. 57] We recall the system of Example 9 (which is equivalent to (12)) with 
iA = s: 


[^J 


d 


dF = 


17 


0 10 
0 0 1 
£4 0 X 


F. 



































As mentioned before, we have a turning point at x = 0. A () 0 is nilpotent and the eigenvalues of 
Aq(x) are 0 and x, whence s = 1. Let T = diag(l, x, x 2 ) then F = TG yields: 


d 



0 

1 

o' 


0 

0 

0 


0 

0 

0 

II 

o 

0 

0 

1 

+ 

0 

0 

0 

X + 

0 

-1 

0 


0 

0 

1 


1 

0 

0 


0 

0 

-2 


Setting fg—x 3 e, the former can be rewritten equivalently as: 


[^1 


^ f)G 


0 

1 

0 


0 

0 

0 


0 

0 

0 

0 

0 

1 

+ 

0 

0 

0 

ft + 

0 

-x 4 

0 

0 

0 

1 


1 

0 

0 


0 

0 

-2x 4 


The leading constant matrix Ao.o is no longer nilpotent. Hence the system can be decoupled into 
two subsystems upon setting G = TW where 



1 

0 

1 


-1 -1 

0 

T = 

0 

1 

1 

+ 

-1 -1 

-2 


0 

0 

1 


-1 -1 

0 


ft+ 0 (d)- 


The resulting equivalent system then consists of the two decoupled lower dimension systems 
where W = [W\ W 2 ] J , cr B — o~c = cr^ = -3: 


[QJ 


f 2 B x 5 dW 1 = 


fix 5 dw z = {1 + fc + (1 + 2x*)fl + O(fl) } W 2 . 


ft + 


-1 

1 + x 4 

:3\ i tj/2 


fl + 0(fl)}W 1 . 


The exponential part (and hence formal solution) of the second subsystem is clearly J f c 2 x 5 (1 + 
0(fc)) dx = ^e _2 x 2 (l + 0(e _2 x 2 )). One can observe that the eigenvalue x of the leading matrix 
coefficient of the input system [A^] is recovered as expected. This is in accordance with the 
exponential parts obtained for (12). As for the first subsystem, B o(x) and Bqq are simultaneously 
nilpotent. Due to this dual nilpotency, in order to proceed in the formal reduction of [Bug], we 
will make use of the e-rank reduction in the following subsection. 


3.2. e-Rank reduction 

We consider again system [A^] given by (8). We assume without any loss of generality that 
Ao(x) and Ao o are simultaneously nilpotent. In this section, we investigate the s-rank reduction 
of the system, i.e. we seek to determine the minimal integer e-rank among all systems equivalent 
to [Ao-J. If this minimal integer rank turns out to be non-positive then we continue the reduction 
by treating the system as in Section 2.1.1. Otherwise, the minimal integer rank gives an upper 
bound for the e-exponential order, which allows us to proceed to Section 4. 

In analogy to its unperturbed counterpart, we define the e-Moser rank and e-Moser invariant 
of system [A , r( ] to be the following rational numbers respectively: 


m e (A) 

Fe(A) 


max (0, h + 


rank (A 0 (x)) 


) 


and 


n 

min {mJT[A rrA ]) for all T in GL n ('K)\. 
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Definition 4. System [4, r( ] (the matrix A (x, £ A ) respectively) is called e-reducible ifm e (A ) > p e (A). 
Otherwise it is said to be e-irreducible. 

If m e (A) < 1, then h = 0. Hence, we restrict our attention to the case of m e (A ) > 1. 

This definition is not to be mixed neither with the usual sense of Moser-irreducible unper¬ 
turbed system [A]; nor with the sense of reduced system already employed in the literature, i.e. 
the system whose matrix is the leading coefficient matrix of [A^]: x p dF = Aq(x)F. In fact, the 
motivation behind Moser’s work in [32] was to determine the nature of the singularity (regular 
or irregular) of an unperturbed system x p 8F = A(x)F, A(x) e Al„(C[[x]]). Consequently, the 
reduction of the so-called Poincare rank p to its minimal integer value was investigated. How¬ 
ever, given our perturbed system [A (r J, it seems more plausible to reduce li, rather than p, to its 
minimal integer value, since a system with h < 0 can be treated as in Subsection 2.1. 

Since p t: (A) cannot be computed from the outset, we aim in this section to generalize Moser’s 
criterion to detect whether the g-rank of a given perturbed system is minimal among all equiva¬ 
lent systems. We remark that Moser’s criterion has been generalized as well to linear functional 
matrix equations in [ 6 ], and borrowed from the theory of differential systems in [24] to investi¬ 
gate efficient algorithmic resolution of the perturbed algebraic eigenvalue-eigenvector problem. 
However, despite their utility and efficiency for such univariate systems, algorithms based on this 
criterion are not considered so far over bivariate fields. 

In [2], we generalize the work of [9, 32] for the case cr A = 0. Herein, we establish the 
following theorem without any restriction on cr A (see [28, Section 5.5, p. 101] and Appendix 
C for the proof). This yields an algorithm which takes an input system [A CTa ], and outputs a 
transformation T e GL n ('K) and an equivalent system [A, r( ] := '/[A fr J which has a minimal 
e-rank among all systems equivalent to [A aA ]. 

Theorem 2. Consider system ( 8 ) given by 


[A(rJ &x P dF = A(x, UW = Yj F 

k =o 


and suppose that h > 1 and m e (A) > 1. The polynomial 

8 a (A) := &"»* <*>(*» d ettf / + ^ + A^x))^ (15) 

vanishes identically in A, if and only if there exists a transformation F = TG, T e GL n (fK), 
such that the equivalent system [A^] := ^[A^] has either a strictly lower e-rank or a leading 
coefficient matrix with a strictly lower algebraic rank. Moreover, cr T , cr A > cr where 


- P 4- ^ P< ] 

otherwise 


Example 8 . Consider the system [A^] q A dF = A(x , %a)F with cr A = 0 and 


A(x,£ a ) = 


2 xfl 

0 

0 

—2x 


3x 2 f 4 

A A 

0 

0 


lx f~A 

0 

£ 2 

a a 

0 


(2x+\)fl 

0 

0 

0 
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Our algorithm computes the transformation F — TG, where T = 


, with cr T — 


0 0 0 
0 0 0 
^.000 
0 0 0 1 


cr A = 0 which results in an equivalent e-irreducible system [A^] 



1 

0 

0 

- A 

0 


2x 

2xf A 

3x 2 

2 x + 

Mx,f A ) = 

0 

0 

f 2 - 

0 


0 

~2xf A 

0 

0 


f|<9G = A(x, f A )G where 


and cr A = 0. One can observe that the s-rank is diminished by tw’o. 

Example 9. Consider the first subsystem resulting from the reduction of Example 7: 




f B x 5 dw l = { 


0 1 
0 0 


-1 0 

-1 0 


fn + 


1 -1 
1 -1 + x 4 


fl + O(fl) I W l . 


It is easy to see that the leading coefficient and constant matrices coincide and that 0 B (A) = 1. 
Hence, the s-rank which is equal to 2 is minimal. 


To proceed in the formal reduction of Example 9, one needs to introduce a ramification in 
s (resp. f B ). The computation of the necessary ramification will be discussed in Section 4. To 
establish the results therein, the e-Moser invariant of equivalent scalar equations is needed. We 
thus introduce it in Subsection 3.3. Before proceeding however, we discuss the case It - l which 
is not tackled by Theorem 2. 


3.2.1. The case h — 1 

The limitation of Theorem 2 to h > 1 is of technical nature (see the proof of necessity). In 
this subsection, we treat the case h = 1. We consider again the system given by (8): 

co 

[A<rJ f A x p dF = A(x, f A )F = A k (x)f k F. 

k =0 

Let h true denote the minimal integer e-rank which can be attained upon applying a transformation 
in GL n (fK) and <o, : (A).\ denote the e-exponential order of [A^]. In this subsection, we give a 
method to decide whether h true = 1 or h true = 0. We first apply the substitution e = e " +1 to 
[A, r J. We denote its e-rank by h. Since h — n + 1 > 1, we can then apply the e-rank reduction 
of Theorem 2 and we denote by h true and wg the true e-rank of the resulting system and its 
e-exponential order respectively. Under these notations, we have: 

Lemma 3. If htrue = 1 then h true — 0. 

Proof. Let a> B = j where ( and d are coprime natural numbers. Then, on the one hand, d, t > 0, 
d < n and to E < h true . On the other hand, oj b = (n + l)a> e and < h true (the minimal e-rank 
bounds the e-exponential order, see Corollary 2). Hence, we have, 

£ £ l 

(o B — (jt + 1 )u) £ = (n + 1 )— > (n + 1 ) — ^ (1 + —^£. 

d n n 
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Hence, if l + 0 then 


1 < (1 + ) l 5: Ct>£ < htrue- 

n 

□ 

Lemma 3 moves the problem from the case li - 1 to the case h > 1. It can be restated as 
follows: Within the formal reduction, whenever the case h = 1 is encountered, a ramification 
e = e " +1 is applied so that the reduction of Theorem 2 can be applied. After this reduction, one 
of the following two cases arises: 

• If the e-rank of the resulting system is less than or equal to one then h, rue of the original 
system is zero, and so we stop our reduction (the e-exponential part have been computed 
completely). If one wishes to continue reduction, then the classical Arnold-Wasow ap¬ 
proach can be employed: First, put the system in Arnold’s form [3] and then use the 
Iwano-Sibuya’s polygon to determine the proper polynomial transformation required to 
arrive at an equivalent system whose e-rank is zero [23], The reduction then continues as 
explained in Subsection 2.1.1. 

• Otherwise, we proceed with the reduction: If the leading coefficient matrix of the ramified 
system has two distinct eigenvalues, then we treat turning points (if any) and apply splitting 
lemma. Otherwise, we proceed to Section 4. 

We remark that in the implementation, we first try to find whether there exist constant vectors in 
the left null space of Ga(A ) (see Appendix C). If such vectors do exist, we use them to construct 
a transformation which might reduce the rank. If the e-rank remains one after all the constant 
vectors are exhausted, we compute a candidate for the ramification in £4 from the characteristic 
polynomial. If this candidate, after applying again the e-rank reduction results in a system with 
non-nilpotent leading coefficient, we keep it. Otherwise, we use Arnold-Wasow approach. 

3.3. e-Moser invariant of a scalar equation 

Due to the equivalence between a scalar differential equation [a an ] and its companion system 
[A (r f\, it is natural to define and question the e-Moser invariant of the former, in the hope to 
gain more insight into the problem. This is the goal of this subsection which is fulfilled by 
generalizing the analogous notion discussed for unperturbed scalar linear differential equations 
in [32, Part IV], 

We consider again the singularly-perturbed linear differential equations given by (9) 

[<V„] d n f + a„-i(x,f a ) d n ~ l f + ... + afx,f a )df + a 0 (x, £,)/ = 0, 

where a n (x,f a ) = 1. We prove the following proposition: 

Proposition 3. Given the differential equation [a (r f. Let r, v be the smallest integers such that 
val E ( afx,f a )) > (i - ri)(j — 1) - v. In other words, let 

k = min [p € N| val e (a,-) + (n — i)g >0, 0 < i < n) 

v — max {(/ - h)(t - 1 ) - val e ( a ,), 0 < i < n) 

Then the e-Moser invariant of the system given by the associated companion matrix is 

v 

p B (a) = *■+-. 

n 
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(16) 

(17) 



Figure 4: Geometric Interpretation 


Remark 2 (Geometric interpretation). Consider again the e-polygon N E (a ) o/[<VJ constructed 
in Subsection 2.1.2) in a ( U , V)-plane. First, one can construct the straight line passing through 
the point (n, 0), with the smallest integer slope k, which stays below N B (a). It is the straight 
line of equation V = k(U - n). Then, one finds among all the parallel lines of slope (k — 1), 
the highest straight line with an integer V-intercept which stays below N e (a). The V-intercept 
of the latter is —v — ii(k — 1) where v is an integer. In other words, the latter has the equation 
V = (U - h)(k - 1) - v. 


Example 10. Let cr a - — 1 and if , — xF a e. We consider the scalar differential equation 


[^tr n ] 


d 5 f + (xf~ 3 + x)d 4 f + 2 x 2 f~ l d 3 f + (f; 3 + 


1 )d 2 f + (~3fa 4 + X 2 f~ 2 )df - f- 2 f - 0. 


First, we plot the points (i, vaff af), 0 < i < n — 5 and the e-polygon (see Figure 4). Next, we 
plot the straight line passing through the point (n, 0) = (5,0) with smallest integer slope such 
that it stays below all these point. This yields V = 3 U - 15 = 3(U — 5) and t = 3. Then, we plot 
the straight line of slope t — 1=2 which stays below these points: V = 2 U - 11 = 2(U — 5) — 1. 
This yields v — 1. Finally, we consider y — ( yo ,.... yf) = (-11, -9, -7, -5, -3). We observe that 
j'o = n — v = 4 and so the equivalent system is given by: 


IArJ 


xf\dW = 


nfi 

0 

0 

0 


X^A 

9^1 

0 

0 


3 £7 


3x? a - x s f 


0 

X^A 

7 & 

0 


0 

0 

xfA 

5 e A 

-2 


0 

0 

0 

xfA 

-x 2 - (x 2 - 3 


W, 


with (t\ — cr a — — 1. One can verify that the system [Ao- a ] is e-irreducible since 0 a(T) = xA 4 . 
Proof. (Proposition 3) For clarlity, we set f := if ,. We define 


y/ = max {x(i - n), (k - 1 )(i — n) — v|, 0 < i < n — 1 . 
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Then, by construction, y,- < val E (a,) and equality is attained for at least one i>n-v (otherwise 
v and k can be minimized). Let represent the smallest integer such that y io = K(i tt - n), then 
j'o = n - v. Geometrically, the y,’s represent a broken line dominated by the (i, val E (a,)). They 
will aid in the construction of an e-irreducible system which is equivalent to (9). In fact, let 

Wi +1 = d'f, 0 < i < n - 1. 


then we have 


dw i+ i 


Uw i+2 

\w i+2 


0 < i < k) 

iq < i < n — 2 


+ 

^ X 


W,+ |], 


‘ W„ - d n - x f - f ‘ 3"- 1 /, 

^ w„ + f ^ + Z-Jo <*i(x,s)w i+ 1]. 

where <r,(x, £) = -a,(x, «, ,k(x)t; k ~ 7i . Let W = (wi,... w„) r and A(x, £) = 


o- a y 0 ^ x£ 


CaTio-l? 

c a yi 0 t K x 

O-aTh+l? 

X 

xa 0 xai ... xa io xa io+ i ... xa„_i + cr a y w _i£* 


Then one can verify that, 

[Aa-J xf K dW = A(x,f)W. 

It remains to prove that this system is e- irreducible. We first remark that A(x, 0) has rank v = n-io 
due to the linear independence of its last v rows. In fact, it’s clear that it has v - 1 linearly 
independent rows (the rows with the x’s). Moreover, a,(x ,£) = -a,(x, £)£~ r ' = -a,{x,f)g~ vals ( ai \ 
and so ai(x, 0) + 0 for at least one i e {0,..., L). Thus the last v rows are linearly independent. 
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Setting f = x‘ T " s we have: 


9a(A) - C det(/i/ H-—— + Ai(x))|^ = o 


A + 6(cr a jo) x 


A + 6(cr a y, 0 _i) 


d + 6(o- a y io ) 


/i + b(cr fl y, 0+1 ) 


-yo'o 

£ 


XQ 1 ! 

£ 


A + <5(o- a y„-i) + 


gv-n+iQ 


A + 6(cr a yo) x 


A + 6(cr a y, 0 _,) x 


xa o 


xa i 


xai, 


where 6 — 1 if k = 1 and is zero otherwise. Hence, since a;(x. 0) + 0 for at least one i e 
{0,..., /o}, 6a(A) does not vanish identically in A and its highest possible degree is iq. It follows 
from Theorem 2 that the system is e-irreducible and p E (a) = p r .(A) - k + □ 

Corollary 1. Given an e-irreducible system [A, Ea ] f h A x p 8F = A(x,La)F. Let [a (Ta ] be the 
scalar equation corresponding to a companion system equivalent to [A^]. Then for k and v 
defined in (16) and (17) we have: k — h and v = rank(Ao). 

Proof By the e-irreducibility of [A^] and its equivalence to the companion system, we have: 

v r 

K + ~ = Ps(a) = p e (A) = m E (A ) = h + - 


Consequently, k = h and v = rank(Ao). 

Corollary 2. Under the notations of Corollary 1, we have: 


v rank(Ao) 

k - 1 h— < a) e {A ) < k and h — 1 -I -< a> e (A) < h. 

n n 


Proof. Follows from (16), (17), and (13). 


□ 


□ 


4. Computing the e-exponential order 

We consider again system [A^] which is given by (8): 


[Aa J x>f h A dF = A(x^a)F = YjAkix)^ F. (18) 

k =0 
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In this section, we assume without loss of generality that h > 0, [A, r J is e-irreducible, and 
,4 (l (x), /4 ( )o are both nilpotent. The goal of this section is to compute the e-exponential order 
to l: (A) and e-polynomials of the system, which give indispensable information within formal 
reduction. In particular, the former determines the ramification in e which can lead to a system 
whose leading matrix is non-nilpotent, so that the process of formal reduction can be resumed. 
This leads to the recursive Algorithm 1. Eventually, we can construct fundamental matrices of 
formal solutions in any given subdomain. We illustrate this algorithm by examples and motivate 
items for further investigation in Appendix E. Let 

A (x s') 

det(T7- "-jr-) = A" + a„-i(x,( A )A„-i + ■ ■ ■ + a 0 (x,f A )- (19) 

XP $A 


such that a„ - 1 and a,(x,f A ) = aij(x)f' A 6 X for 1 6 {0,We define the e- 

polygon N e {A) of [4 (TA ] by taking P e (A) of Section 2.1.2 to be the union of P(i, val e (ai(x,^ A )) 
for i e {0,..., n). We can then prove the following: 

Theorem 3. Consider the s-irreducible system [ A a - A ] given by (18) with h > 0 and (19). If 
h > n — rank(Ao(x )) then the e-exponential order o/[Ao- a ] is given by 


w e (A) = max ( 

0 <i<n 


-val e (aj ) 
n — i 


and its corresponding s-polynomial is given by the algebraic equation 


e 

E e (X) = y]^ AVaUaik) a ik ,vaid*, t ) x(it ~ io) 
k=0 

where 0 < io < i\ < ■ ■ ■ < ig = n denote the integers ifor which a> E (n - i) = —val £ (aj) (i.e. lie on 
the edge of slope cj e of the e-polygon N e (A ) o/[4„.JJ; and a ai(x,tj A )\ u=0 . 

The proof can be established by adapting to the parametrized setting the proofs of [ 8 , Lemma 
3, Lemma 4, Proposition 1, Theorem 1] (see Appendix D). Not only does this theorem compute 
these invariants of the system, but it also allows a further reduction of the system as follows: Sup¬ 
pose that a> E (A ) = ^ where f d are relatively prime positive integers. One can then set e = e 1 ' d 
(or equivalently 3 f A = x (r(d ~ 1 >ld ) in [4 cr J and perform the e-rank reduction (the minimal e- 
rank is I). This will lead to an equivalent system whose e-rank is f and whose leading coefficient 
matrix has at least d distinct eigenvalues. The system can thus be decoupled into subsystems of 
lower dimensions. By repeating these procedures for each of the resulting subsystems, we can 
decouple the initial system into subsystem(s) of lower dimension(s) or zero e-rank. This leads to 
the recursive algorithm of Section 5. 

We remark that the condition h > n — rank(Af x)) in Theorem 3 is non-restrictive. We can always 
arrive to a system satisfying this condition by a generalization of [ 8 , Lemma 5], which is based 
on applying the affinity (U,V) —> (U,d V ) for some integer d to the e-polygon. In fact, we can 
prove the following lemma: 


3 In fact, we want to set eA = e. And so we choose an integer s such that f\ ‘ 
yrrfgrf = x°- +s e. 


£ax s which is equivalent to 
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Lemma 4. Consider an s-irreducible system [A^] given by (18) with h > Oand r — rank (Ao(x)). 
Let d be an integer such that d > A _” +£ and let 

i A h x p dG = A(x, £a)G. 

be the s-irreducible differential system obtained by the ramification s = sf 1 (or equivalently 
£a - f A d x crA ^ 1 ^ d and cr A — cr A (d - 1)) and performing s-rank reduction. Then pffA) = h+ L , 
where r = rank (Af x)) and h + r > n. 

Proof. By Corollary 2 we have, h — 1 + ^ < 0 J r (A) < h. And due to s — sf 1 , we have u> E (A) = ‘AdAl. 
Hence, w g (A) > d(h - 1 + jj) > n. But h > ntg(A) and r > 1, which yields h + r > njg(A) + 1 > 
d(h - 1 + '-) + 1 > n + 1 > n. □ 

5. Formal reduction algorithm 

With the algorithms of Splitting lemma, turning point resolution, e-rank reduction, and as can 
be verified by Algorithm 1 below, we have re-established constructively the well-known general 
form for a fundamental matrix of formal outer solutions for an input system [A £ ] given by (3): 

oo 

F = (2 <f>*(x 1/s ) (xf s) kld ) exp( 
k=0 

where s, d are positive integers; cr is a nonpositive rational number; <3 is a diagonal matrix whose 
entries are polynomials in s l/d with coefficients in C((x'^ s )). We refer to Q := f Q as the s- 
exponential part (logarithms in a root of v might arise as a result of integration); and the entries 
of the 1 ^'(’s are root-meromorphic in x (see [36, Introduction] or [35]). 

Remark 3. Under the notations and statements of Theorem 3, the leading term of Q is given by 

J diag(Ai (x),..., X deg (lE Jx ), 0,..., 0) dx, 

where the Xfs denote the roots of the s-polynomial E e (X). A similar statement can be stated in 
terms of the eigenvalues of the leading matrix coefficient, x p , and f h A . 

We sum up our main results in the formal reduction algorithm. Algorithm 1, which computes 
the s-exponential part and consequently a fundamental matrix of formal solutions (20) in a given 
subdomain. We recall that: 

• If n = 1 then we proceed by integration up to the first h terms; 

• If /z < 0 then we follow Subsection 2.1.1 (using ISOLDE, miniIsolde, or Lindalg). 

6. Conclusion and further investigations 

In this article, we give an algorithm, implemented in Maple 4 , which computes a fundamental 
matrix of outer formal solutions for singularly-perturbed linear differential systems in a neigh¬ 
borhood of a turning point. The subprocedures discussed are stand-alone algorithms (Splitting, 


/ 


<3(x 


1 Is _-l Id 


)), 


( 20 ) 


4 The package is available at: http : / /www.specfun.inria.fr / smaddah l Research Jit ml 
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Algorithm 1 Exp_param(/z,p, ct a ,A(x, £ a )) : Computes the e-exponential part of system [A, rA ] 
(performs formal reduction) 

. Input: h,p,cr A , A(x, £ 4 ) of system ( 8 ) (p - cr a - 0 for system (3)) 

Output: Q(x^\s~ ]/d ) 

Q <- diag(0,... 0); 
while h > 0 and n 4- 1 do 

if Ao q has at least two distinct eigenvalues then 
Apply the e-block diagonalization of Section 2.3; 

Exp_param(/i, p, (t a ,A 11 (x, e)); Update Q and <x /4 ; 

Exp_param(/z, p, ct a ,A22(x, e)); Update Q and cr A \ 
else if Ao.o has one non-zero eigenvalue then 
Update Q from the eigenvalues of Ao.o; 

A(x, £ 4 ) <— perform eigenvalue shifting (14); (Ao.o is now nilpotent); 

Exp_param(/i, p, <t a ,A(x, £ 4 )); Update Q ; 
else if Ao(a') is not nilpotent then 

A(x, £ 4 ), p, cr A , <— apply turning point resolution of Subsection 3.1; (Ao.o is now non- 
nilpotent and cr A is updated); 

Exp_param(/i, p, ct a ,A(x, £4)); Update Q ; 
else 

A(x, <? 4 ), h <— e-rank reduction of Section 3.2; 
if /: > 0 (i.e. h true > 0 ) then 

if Ao.o has at least two distinct eigenvalues then 
Apply the e-block diagonalization of Section 2.3; 

Exp_param(/;, p, <r A . A 11 (x, e)); Update Q and cr A , 

Exp_param(/z, p, cr A , A22 (x, e)); Update Q and <x A ; 
else if Ao.o has one non-zero eigenvalue then 
Update Q from the eigenvalues of Ao.o; 

A(x,£ a ) <— perform eigenvalue shifting (14); (Ao.o is now nilpotent); 
Exp_param(/z, p, ct a ,A(x, £ a )); Update Q \ 
else if A()(x) is not nilpotent then 

A(x, £ 4 ), p,cr A ,<— apply turning point resolution of Subsection 3.1; (Ao.o is now 
non-nilpotent and cr A is updated); 

Update Q from the eigenvalues of Ao.o; 

A(x,£ a ) <— perform Eigenvalue shifting; (Ao.o is now nilpotent); 

Exp_param(/z, p, ct a ,A(x, £ a )); Update Q\ 
else 

Use Theorem 3 of Section 3; 
co E - |; £ <- E d \ 

/i, A(x,£ a ) <— Apply e-rank reduction of Section 3.2 (h <— C)\ 

Exp_param(/z, p, ct a ,A(x, £4)); Update Q; 

end if 
end if 
end if 
end while 
return (Q). 


27 





rank reduction, and resolution of turning points). The formal reduction algorithm is based on the 
generalization of the algorithm given in [8] for the unperturbed counterparts of such systems. 

Our results are presented in the formal setting. However, the growth of the order of poles in x 
is tracked within the formal reduction. This gives information on the restraining index of the sys¬ 
tem, and consequently, an adequate stretching transformation can be chosen and performed. This 
furnishes the first step in resolving Iwano’s first problem and computing, for an input system, [p] 
which is given by (11): 


[p] 0 = po < pi < p 2 < • • • < p m - 

In Appendix E, we try to motivate the employment of our proposed formal reduction algo¬ 
rithm to resolve the same problem for a general system: Suppose that we start the reduction with 
the input system [Aq-J, <. r ,\ = 0, given by ( 8 ) as an input. Then, upon applying Algorithm 1, we 
can determine the outer solutions and the final restraining index cry,„ fl/ which allows full reduc¬ 
tion. If the restraining index is nonzero then we have a turning point and p m = -l/o-f ina i. We 
can then perform the stretching r = xs~ Pl in the input system [A cr J, and apply to it Algorithm 1. 
This determines p 2 - We show an example where the same process can be repeated with p 2 to 
determine p 3 , and iteratively we reach a system whose final restraining index is zero. In future 
work, we hope to investigate this approach and its correctness. 

Another field of investigation is the relation between the algebraic eigenvalues of the matrix 
A(x, s ) of such systems and the exponential part of the solution (see Example 21 and [28, Chap¬ 
ter 2]). One can then benefit from the existing work on fractional power series expansions of 
solutions of bivariate algebraic equations, to compute the restraining index (see, e.g. [17, 4] and 
references therein). On the other hand, it would be interesting to investigate a differential-like 
reduction for a two-parameter perturbation of a JCF (see [28, Chapter 2] and references therein 
for the one-parameter case). In fact, one can observe that the main role in the reduction pro¬ 
cess is reserved to the similarity term of 7’[A 0 - /( ], i.e. T 1 ,4 rri 7' rather than T l dT. Hence, for 
a non-differential operator where T [A lfA ] = T 'A tr J\ the discussion is not expected to deviate 
substantially from the discussion presented here for a differential one. 

The generalization of other notions and efficient algorithms can be investigated as well (see, 
e.g. simple system [13, 11] and [28, Appendix A]). 

In Examples 18 and 20 we treat input systems with cr A < 0. The algorithms generalized 
herein have been generalized to treat systems with an essential singularity [ 12 ] and difference 
systems as well(see, e.g. [1, 6 ]). This motivates investigating the adaptation of our proposed 
algorithms at least to the difference setting. 

There remain as well the questions which fall under the complexity, e.g. studying the com¬ 
plexity of this algorithm; the bounds on number of coefficients needed in computations; devel¬ 
oping efficient algorithms for computing a cyclic vector; and comparing it to the results of our 
direct algorithm. 

And finally, to give a full answer on the behavior of the solutions, the related problems of 
connection, matching, and secondary turning points, are to be studied as well (see, e.g. [ 20 , 22 ] 
and references therein). 
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Appendix A. Computing a companion block diagonal form 

It is well-known that a scalar equation [fvj given by (9) can be expressed as a first-order 
system by setting 

F = [f,df,d 2 f,...,dr- l n T . 

However, the opposite direction of this transition is nontrivial although possible. The theoret¬ 
ical possibility stems from the work of [19] with the so called cyclic vectors. Consider again 
system ( 8 ) given by: 

[Ao-J $ h A x PA dF = A(x,% a )F, = x (TA e and A 0 , 0 * 0. 

We also recall that with the definition of £ 4 , A(x,£a) = Yik=o A*(jc)^ with A^.(x) e A1„(<C[[x]]) 
for all k > 0. In this appendix, we compute a companion block diagonal form for [A fr J, in a finite 
number of steps, in a neighborhood of a turning point. This algorithm is an adaptation of that 
of [7], developed for the unperturbed counterpart of [A, r J. It relies on a sequence of polynomial 
(shearing) transformations in x, and elementary row/column operations. 

• Shearing transformations: 

T = diagC^ 1 ,* 0 - 2 ^ 2 ,.. .,**•#) (A.l) 

where the cr’s and /j’s are respectively rational numbers and integers. We remark that 
shearing transformations may alter cta because of the a’s. 

• Elementary transformations: We consider transformations of the form 

T = I„ + Pij(a) (A.2) 

where P, , is a zero matrix except for entry in the i th row and j th column for i, j e [ 1 ,..., 11 } 
and a e R such that a(x = 0, £4 = 0) 4- 0. Obviously, T € GL„CR). The equivalent system 
[Ao-,] = TtAo-J is such that h = h, p = p, and cr^ > min(^ + cr A , cr A ). Moreover, it can be 
easily verified that the effect of this transformaiton on the system [A, r J is as follows: 
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Mfa ): Multiplies the i' h column and i ,h row by a and 1/a respectively. It also adds 
to the diagonal entry in the (i, i) th position. 

Cjj(a): Adds to the j ,h column the i' h column multiplied by a, adds to the i th row the j' h 
row multiplied by -a, and adds (—g A x p ida) to the entry in the (i, j) th position. 

Rjj(-a ): Adds to the i ,h row the j th row multiplied by a, adds to the j th column the i th 
column multiplied by -a, and adds {f h A x p da) to the entry in the (i, j) th position. 

A careful choice of these transformations will allow us to establish the following result construc¬ 
tively: 

Theorem 4. Consider the system |A (r J given by ( 8 ). Then there exists a transformation T 
which is a product of transformations of the forms (A.l) and (A.2) in a root of x, such that the 
equivalent system [A^] ^Ix^dF = AF has the block diagonal form: 

A(x,f A ) = diag (A u (x,f A ),... ,A {e (x,f A )) 

where the An,..., A cc are companion matrices. 

The proof of this theorem relies on two lemmas which we establish herein. Before pro¬ 
ceeding, we can always assume, without loss of generality, that Ao.o = A (x = 0, £ A = 0) is in 
Frobenius canonical form, i.e., 

Ao,o = diag(A^ 0 ,...,A«) (A.3) 

where the block submatrices A^,... ,A e Q e Q are constant matrices in companion form, of dimen¬ 
sions n\ > n 2 > ■ ■ ■ > H{ > \ («i + «2 + ■ ■ ■ + >i{ — n). 


Lemma 5. Consider the system [A^J given by ( 8 ) whose leading constant coefficient Ao.o is in 
Frobenius canonical form. Then, by a sequence of elementary operations of the form (A.2), the 
equivalent system [A^] f l lx p dF = AF has the following block decomposition: 


Mx.f 4 ) = 


S(x,f A ) 

V(xffi A ) 


U(x,f A ) 

W(x,f A ) 


(A.4) 


where 


• S (x,f A ) is a n\-square matrix in companion form: 


S(x,t A ) = 


0 

0 


1 

0 


At $2 


0 


1 

Sn\ 


• U(x,^ 4), V{x,^f) are of dimensions n\ X {n — n\) and (n — n i) X n\ respectively, Uo,o = 
Onix(n-n\)> Vo,o ~ O (n-n\)xn\ > cind of the forms. 



0 . 

0 


v\ 0 . 

. 0 

U(x,f A ) = 

0 . 

U\ 

0 

• U n -n\. 

and V(x,t A ) = 

yn-ti\ 0 

. 0 
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W(x, if A ) is a square matrix of dimension n - n \; 


• Furthermore, Aq.o = Ao.o- 


Proof. We first assume n\ > 1. Otherwise, A(x, if A ) = [a;/] 1 <;,/<« is already in form (A.4). Since 

Aoo has the form (A. 3), a,j +1 (x = 0,f - 0) = 1. Hence, 1 /a l>i+ i (x, if A ) e % and 1 /«,■;+1 (x = 

0,fr =0)= 1. 

Step 1 For each row number i from 1 to n\ - 1: 

(1.1) We use l/ai,i+i(Jt,£i) as a pivot and apply M, + i(l/a Ijl+ i(x, £ 4 )) to set the entries in 
the (i, i + 1) positions to 1. We remark that the term (-iffx p (a)d( 1 /a)) is added to the 
diagonal entry in the (i + 1 ,i+ l)' h position. 

(1.2) We now make use of these 1 ’s to set the entries in the (i, j ) positions for i e {1,..., n\ — 
1), j e {1,..., n}, j + i + 1 to zeroes: For each row number j 40+1 from 1 to n, 
we apply C,+i j(-a,j). We remark that the term (+if h A x p da) is added to the entry in 
the (i + 1, j) position. Clearly, the rows 1 to f — 1 are not altered by these elementary 
operations. 


Step 2 We have attained the anticipated form for the first n ] rows, i.e. we have shown how 
to construct an equivalent system whose first n 1 rows have the described properties of 
[S(jr,£ 4 ) U{x,f A )\. If n j < n then we can proceed to work on the last n - n\ rows. 


(2.1) Again, we make use of these l’s created in Step [(1.1)] to set the entries in the (i, j ) 
positions for i e {n 1 + 1and j e [2,... ,n\) to zeroes: For each i, j in these 
ranges, we apply successively for each i from n\ + 1 to n with an inner 

loop of j from ti\ to 2. Evidently, each such operation does not alter the form created 
by Step 1 for the first n\ rows, and by its preceeding operations within Step 2. 

Finally, since each of the operations performed is of the form (A.2), for x = 0 and £4 = 0, such 
an operation reduces to Hence, the form of ,4 0.0 is preserved. □ 


Remark 4. One can observe that by a repetitve application of Lemma 5 on the diagonal blocks 
of dimensions ti\ > n 2 > • • • > ti(, we arrive at the aformentioned Arnold form. 

Lemma 6 . Consider the system [Agiven by (8). Suppose that its leading constant coefficient 
Afyj is in Frobenius canonical form and A(x, iff) is inform (A.4). Then, by a sequence of constant 
transformations, elementary operations of the form (A.2), and shearings of the form (A.l), the 
equivalent system [A^] ^Ix^dF = AF is such that has the following block diagonal form: 


Mx,fA) 


S(x,f A ) O 
O W(x,f A ) 


(A.5) 


the inner loop 
shld start from 
because the 
term coming 
from the deriva¬ 
tion settles in 
the (i,j - 1 ) 
position and so 
will be deleted 
afterwards. I 
should specify 
how cr A changes 
w.r.t. cr A . 


where S(x,f A ) is a r-square companion matrix where r > n\ and W(x,f A ) is an - r square 
matrix. 


Proof. If IJ(x, La) and V(x, if A ) of (A.4) are smiltaneously zero matrices or n \ = n then A(x, 44) 
is already in the form (A.5). Otherwise we suppose that U(x,f 4 ) is nonzero and limit our 
discussion to this case. If V(x, if a) is nonzero as well, then we can proceed analogously. Let 
uj = up< ( {x)if k A for j e [1 ,..., n - n\ } then we define 


/3 = mini <j<n~ n f val e (uj)) and a - mmi<j< n . n fval x (ujp)). 
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We remark that /? > 0, and it’/) = 0 then a > 0 since £/o,o is a zero matrix. Let 


T = diag(x“^/ ni , and [Bo-J = TTA^J. 


Then we have: 


«(*,&) = 


S(x,t A ) - xr +al f A + \a+/3o- A )I ni 
x a f A V(x,f A ) 


X- a f/U(x,tA) 

W(x,Z A ) 


and 


Bo(x ) = 


5 0 (x) 
O 


x- a f/u 0 (x) 

W 0 (x) 


B o,o - 


5 0,0 

0 


0,0 

Wo.o 


where S o,o is a companion matrix of dimension «i, Uo.o is a companion matrix whose last row 
has at least one nonzero element, and Wo.o is a block diagonal square matrix whose blocks are in 
companion form. Let e, denote the i' h row of the identity matrix /„ and B' Q () denote the i th power 
of fi 0 ,o. W then have: 


e iB' 00 - ^ Cijej + e,+ \ for i = 1,... ,m - 1 

j= i 

n\ n-n\ 

^o.'o — ^ c ni,j e j + ^ uj,[3e ni+ j 

i =i j =i 


where Cy 6 C. Since there exists at least one j e {l,...,n-«i} such that u^p £ 0, then the 
vectors e\,e\ Boo, e\ B^ 0 , e\ B”' 0 are linearly independent. Hence, denoting by r the dimension 
of the first block of the Frobenius form of Bo,o, i.e. the degree of the minimal polynomial of 
B 0 ,o, we have r > n\. Hence, by putting the leading constant matrix in Frobenius form, applying 
Lemma (5), and the procedure described herein, we either arrive at the desired form in a finite 
number of steps. □ 

Proof. (Theorem 4) The proof can be attained by a recursive application of Lemmas 5 and 6. □ 

Remark 5. We remark that this method can give a cyclic vector and that [A^] is not uniquely 
determined by [A^]. 


Appendix B. Proofs for Section 2.3 

Proof. (Proposition 1) Let f(x,f) = exp(f q(x, f)dx) with q{x, f) e (Jjsn* < C((x))((e 1/i )). We 
remark first that 

Vie N, d'f — Pj(q)f 

where 

Pfq) = <?' + ••• + d‘f l q 
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is a polynomial in (< q , dq ,..., d l l q). Then, /(*,£) verifies (9) if and only if q(x,£) is a solution 
of 

n 

^ai(x y g)Pi(q). 

1=0 

We denote by e the polar order of q in £ (e-order), and write 


q(x, g) = g e (s(x) + 0 (g)), q t 0 . 


Then, Vi e {0,. val E (Pf) - -ie and in fact Pfq(x, g)) = g ,d (s(x)‘ + 0(g)). 

Thus, 

n n 

y^a^gW.iq) = 

i= 0 i=0 

where val E (Wi ) > 0. Let v(e) = mmo<,<„(v, - ie) and 7(e) = {/ G {0, ..., ii) | val e (aj - ie) = v(e)}. 
Hence, 

n 

X ai(x,g)Pi(q) - r W (X «/,v,(x)V W + 0(0). 

!=0 ie/(e) 

Thus, q(x,g) = g e (s(x) + 0(g)), q i 0 satisfies £"=o afx, g)P fq) = 0 if and only if s(x) is a non 
zero solution of the equation 

X a i.v,(x)X‘(x) = 0. 

ie/(e) 

This equation has non trivial solutions if and only if \I(e)\ > 1, i.e. if and only if e e (ei ,..., e/j. 
And then, A(x) is one of the nonzero roots of the associated determinant equation. □ 


Proof. (Theorem 1) We proceed in two steps: 

Step 1: We first block-diagonalize the leading matrix coefficient: Suppose that there exists 
Tq(x) g G/„(C[[x]]) s.t. the substitution F = Tq(x)G in [A^] yields an equivalent system 
[Ao-J g h A x p dG = A(x,g A )G where -.A(x,g A ) = 2^=0 Mx)g k A , A k (x) 6 M„(C[[x]]) for all k > 0, 
and the leading coefficient matrix Ao(x) is non-zero and has the following block-form in accor¬ 
dance with Aq,o: 


A 0 (x) = 


a ( ;'w 

O 


O 

Af(x) 


Then we have: 

A T 0 (x) - T 0 (x) A = g h A x p dT 0 (x) 


which yields 


A k (x) Tq(x) - Tq(x) A k (x) = 6ki,dTo(x), where Ski, = 1 if k = h and Ski, = 0 otherwise. (B.l) 


In particular (since h > 0): 


A 0 (x) T 0 (x) - T 0 (x) A 0 (x) = 0. (B.2) 

We further assume that Tq(x) e A1„(C[[x]]) is in the following form 

I A 2 (x) 

T 2l (x) I ■ 
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T 0 (x) = 






Then, for 1 < g £ g < 2, we have: 


Af(x) - Af(x) + Af(x) Tfix) = 0 
Af(x) + Af(x) Tf(x) - Tf(x)Af(x) = 0 


Inserting the series expansions Ao(x) = 2” 0 A,;ox' and Ao(x) = in (B.3), and equating 

the power-like coefficients, we get for i = 0: 


4® _ see 
l^o.o ~ ^0,0 

) T gg — T® 4 ee — (1 
l^O.O 1 0,0 1 0,0 ^0,0 _ w 


which are satisfied by setting Ail = ASS and TfL = /, Tf' = O. And for i > 1 we get: 


i,0 


i,0 


/-I 

- V(A^. n r^-r^A^. n ) 

Z_J V <~J, 0 ], o J. 0 i-y,0' 

7=1 

i-1 

(B.4) 

V A eg T ge 

Zj i-j, 0 1 /0" 

7=1 

(B.5) 


It’s clear that (B.4) is a set of Sylvester matrix equations that possess a unique solution due to 
the assumption on the disjoint spectra of Ail and A 22 () (see, e.g., [5, Appendix A.l, p. 212-213]). 
Remarking that the right hand side depends solely on the A^o, 7jo, A jk) with j < i, equations 
(B.4) and (B.5) are successively soluble. Hence, such Tq(x) e G/„(C[[x]]) can be constmcted. 
Clearly, pr = v t - 0 and cry- > cr A . Moreover, it follows from (B.l) that v A = 0, p A = 0, and 
cr A > cr A . 

Step 2: Due to the first step, we can assume without loss of generality that the leading coefficient 
matrix Ao(x) of the system [A, Ta | has a block-diagonal form in accordance with that of Ao.o- It 
suffices to seek T and A in the form T = E 4 L 0 ^ki x )i; and A = E 4 I 0 A k (x)% respectively where 
T k (x),A k (x) e GL„(C[[x]]) for all k > 0. We set 7 q(x) = /„, Aq(x) = Aq(x), and we rewrite 


00 00 

A = ^ A k (x)^A = ^ B k (x)Z, where B k (x) = x ((TA_<r) %(x). 
k=0 k=0 


Since cr A > cr, it follows that B k (x) € AI„(C[[x]]) for all k > 0. With the expansions in we get 
the following set of ansatz for k > 0: 


Y(BUx)Tiix) - Tiix)A k _iix)) = x p (d + — — —)T k _ h , where T k _ h = 0 for k < h. (B.6) 

X 


It then follows from (B.6) that for k > 1 we have: 


B 0 (x)T k (x) - T k (x)A Q (x) 


k-1 

= ^ Ih-iixfrj(x) - T\ix)A k -iix) 
i= 1 

„ 

+ x p (d +- )T k -h{x) - A k ix) + B k ix). 
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By inserting the claimed block forms of T k (x) and A k (x) for all k > 1 in the above equation, with 
1 < Q + g < 2, we obtain: 


Af(x) 


B*(x)T*(x) - Tf(x) Af(x) 


Bf(x) 

Y Bfjx) 7f'(x) - ifjx) Af(x) 

i= 1 

X p (d + ^^)7^(x) + Bffx). 


Hence, for every k > 1, the computation of T k (x) requires only lower order terms of T. B, and A. 
Hence, by setting T k (x - 0) - 0 and expanding each of the equations for a fixed k > 1 w.r.t. x, 
we arrive at constant Sylvester matrix equations, which can be resolved successively due to the 
disjoint spectrum of Aqq and A“ 0 , to compute T k (x) up to any desired order in x. □ 


Appendix C. Proofs for Section 3.2: e-rank reduction 

In this section, we establish our constructive proof of the necessity of Theorem 2. For suffi¬ 
ciency, one can consult [28, Section 5.5, p. 101]. We prove the following: 

Theorem 5. Consider the system (8) given by: 

oo 

[Ac -J x P t h A DF = A(x,Z A )F = Yj A k(x)tA k F 

k=0 

with h > 1 and m e (A ) > 1 and set rank (Ao(x)) = r . Suppose that the polynomial 

9a(A) := f A det(/17 H—^ + Ai(x))[f A= o 

vanishes identically in A. Then there exists a transformation F = R(x, £'k)G such that the equiv¬ 
alent system 

OO 

[AaJ t/xPdG = A(jc,ft)G = Y A ^a G 

k=0 

satisfies either h < h or h — h and rank(Ao(x )) < r. Moreover, such a R(x, can be constructed 
by [2, Algorithm 1 ] where it is always chosen to be a product of unimodular transformations in 
GL„(C[[x]]) and polynomial transformations in f A . 

The proof of the theorem will be given after a set of intermediate results. 

Lemma 7. Given the system 

CO 

[Ao- a ] x p f b A 8F = A(x, f A )F = Y ^(4 F, 

k=0 

with r — rank (Ao(x)). There exists a unimodular transformation U(x) in GL„(C[[x]]) such that 
the leading coefficient matrix Ao(x) of the equivalent system (10) given by 

[A ( r A ] f\x p 8G = A(x,f A )G, cr A > min( j- + cr A , cr A ) 
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has the form 


A 0 (x) = 


A“(x) 

Aq(x) 


O 

O 


where A, 1 , 1 ^) is a square matrix of dimension r and 
r. 


A“(x) 

_Al\x) 


is a nx r matrix of full column rank 


Remark 6. In practice, U (x) can be obtained by performing Gaussian elimination on the columns 
of Aq(x) taking as pivots the elements of minimum valuation (order) in x. 

Proof. By Remark 7, Aq(x) = U UxjAfxjUix). Hence it suffices to search a similarity trans¬ 
formation U(x). Since C[[x]] is a principal ideal domain (the ideals of C[[x]] are of the form 
x*C[[x]]), it is well known that one can construct unimodular transformations Q(x), U(x) lying 
in GL„(C[[x]]) such that the matrix Q(x)Aq(x)U( x) has the Smith normal form 

Q(x) A 0 (x) U(x) = diagCx® 1 ,..., A 0,..., 0)) 


where det((/(0)) + 0, det(Q(0)) ^ 0 , and /Si,... ,/8,- in Z with 0 < /Si < /3 2 < • • • < fr¬ 
it follows that we can compute a unimodular matrix U(x) in GL„(C[[x]]) so that its last n - r 
columns form a C[[x]]-basis of ker (Aq(x)). 

As for finding a lower bound of cr A , one can see that since U(x = 0) is invertible, we have: 

, ^ (val x (A k ) \lk + h 

val x (Af) > < 

1 val x (Ai,) + p otherwise. 


Thus, if p < 0 then cr^ > cr A + y r □ 

Remark 7. Consider system [A (Fa ] given by (8). Let T(x) — (/„ + Q(x)) e GL„(C[[x]]) such that 


(J( x) — I q - / ] 1 <i.j<u , 


(qnj G C[[x]] , for r + 1 < j < n 
1 qij = 0 elsewhere. 


One can then verify that the resulting equivalent system (10) is 

\A(r x ] f\x p dG = A(x,f A )G , cr A > o- A 


where 


A 0 (x) = (/„ + Q(x)) 1 A 0 (x) (/„ + Q(x)) 

• Afx) - (/„ + Q(x))-' A 1 (x) (/„ + Q(x)) - yx p dQ(x)\ 
A 2 (x) = ... 


(C.l) 


where y — 1 if h — 1 and y — 0 otherwise. But, due to (15), we can limit our interest to Ao(x), 
Ai(x), Ao(x), and A\(x) solely. Hence, if h > 1 then it suffices to investigate T~ l A(x,£o- a ) T. We 
remark that if h = 1 then the term T~ l dT(x) should be taken into account. This is the reason 
the case h = 1 is treated separately. This distinction is purely technical as explained in the 
following remark and unlike for the Poincare rank p of an unperturbed system, it does not lead 
to a classification of regularity. 
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By Lemma (7), we can suppose without loss of generality that Aq(x) is of the form 


A 0 (x) = 


Ag'(x) 

A 21 (x) 


O 

O 


(C.2) 


with r independent columns and (n - r ) zero columns. We partition A \ (x) in accordance with 


Aq(x), i.e. Ai(x) 


A J 1 (x) A| 2 (x) 
A 21 (x) Af(x) 


and consider 


G a (A) = 


Aq'(x) 

AqCx) 


A| 2 (x) 

A 22 (x) + AI„-r 


(C.3) 


This consideration of G a (A) gives an e-reduction criterion equivalent to 0 a (A). In fact, let DO;) = 
diagl^/,.,/„_ r ) where r = rank (Ao(x)). Then we can write £f A A fx,g A ) = ND~ l where N 
N{x,ff) e M n {K). Set Do = D( 0) and No = jV(x, 0). Then we have 


detlG^U)) 


det(M) + AD 0 ) = det (N + AD) | fA=0 
(det(^^+T/„)det(ZJ))| fA=0 

AA 

(det( 0 + Ai(x) + AI „) ^)|f A= o = &a(A). 


Thus, det(G,i(/l)) = 0 vanishes identically in A if and only if 0 A (A) does. We illustrate our progress 
with the following simple example. 

Example 11. [2, Example 2] Consider [A :rA ] f h A dF — A(x,f)F with cr A = 0, h > 1, and 


A(x,£i) = 


& 

X 2 

—X 

0 


-x 3 £t 

Xf A 

0 

2 


Clearly, Aq(x) is nilpotent of rank 2 and 


G a (A ) = 


0 

x 2 

-x 

0 


We then have the following proposition: 
Proposition 4. Given the system 


(1 +x)£a 
0 
0 
0 


0 

0 

0 

2 


x + 1 
0 
A 
0 


0 

-2xf A 

2& 


0 

—2x 
2 
A 


[A(t a ] ? A x p OF = A(x, f A )F = A k {x)it k F, 


(C.4) 


k=0 


with r = rank (Aofx)). If m E (A) > 1 and det(G, 4 (/l)) = 0 then there exists a finite product of 
triangular matrices T(x) = P(I„ + Q(x)) where P is a permutation and det T(x) — +1, such that 
for the equivalent system 

oo 

[A^] f\x p dG = A(x,&) G = Yj Ak(x) t k A G 

k =0 
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we have: 


GaU) 


Aq 1 (x) U i(x) U 2 (x) 

Vi(x) W\{x) + AI„- r -g W 2 (x) 
_M\(x) M 2 (x) Wfx) + Alp 


(C.5) 


where 0 < Q < n — r, W i, W 3 are square matrices of order (n — r - g) and g respectively, and 


) = rank (^Aq*(x) t/i(x)J), 

rank ([aJ'Cjc) t/i(x)]) < r. 

Moreover, <X| > min(^ + cr A ,cr A ). 

We shall need the following Remark in the proof of the Proposition. 


rank ( 


K'(x) 

Mfx) 


U\(x) 

M 2 (x) 


(C. 6 ) 

(C.7) 


Remark 8. [2, Remark 6] Suppose that Ga (/l) has the form (C.5) and there exists a transforma¬ 
tion T(x) € GL„(C((x))) such that Gt[a„. a ]{A) has the form 


G T[A^ a ]U) 


M 11 

A o 

U 1 

U 2 

Vi 

Wl + AI n - r - 

0 W 2 

0 

O 

W 3 + Alg 


where 0 < g < n — r and W 3 is upper triangular with zero diagonal. Then, 

UV Ui 


detCG^^J^detC^ Wi+ ^ 


)■ 


If det(G 7 -[A ](A)) = 0 then we have det(Gr[A^](d)) = 0 as well (rank of leading coefficient matrix 
is unchanged). Hence, 


det( 


A 11 

A 0 

Vi 


U 1 

W 1 + /\I n — r —Q 


) = 0 . 


(C. 8 ) 


For a fixed g € {0,... ,n — r] we shall denote by G^ the matrix 


4 11 
A o 

Vi 


u 1 

Wi 


of (C.5). 


Proof. (Proposition 4) Since dct(G,i(4)) = 0 then in particular, the matrix Ga(A = 0) is singular. 
Let £1 (respectively If) be the vector space spanned by the first r (resp. last n - r) rows of 
Ga(2 = 0). We have 

dim(E\ + E 2 ) — rank (Ga(A = 0)) < n. 

If dim(E 1 ) < r then setting g = 0 suffices to fulfill our claim. Otherwise, since 


dim{E\ + E 2 ) = dim(E\) + dim(E 2 ) - dim(E\ Pi E 2 ) < n. 


it follows that either dim(E 2 ) < n-r or dim(E \ C\E 2 ) > 0. In both cases, there exists at least a row 
vector m (l \x) = (rn^fx ),..., xzt^\x)) with entries in C((x)) in the left null space of Ga(/ 1 = 0), 
such that m { J\x) 4 0 for some r + 1 < j < n. We can assume without loss of generality that 
vf l> (x) has its entries in C[[x]]. Indeed, this assumption can be guaranteed by a construction 
as in Remark 6 . Let the constant matrix P (1) denote the product of permutation matrices which 
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exchange the rows of A, so that valftal}’ (x)) < val x (v/j\x)) ,r+\<j<n-\, where val x 
denotes the x-adic valuation (order in x ). Let 


Q m {x) = [q] l /(x )], s.t. 


_ rJV 


\&x) = - 


"(x) 


qfi = o 


w 

elsewhere 


for r + 1 < j < n 


Thus, P {X \l„ + <2 (1, (x)) is unimodular in GL„(C[[x]]). 

Set F = F°, A = A (0) , A = A (£>) and let A (1) be the matrix of the equivalent system f h Am x p dF (X) = 
A^\x, ^ 4 (i))F^ 1 - ) obtained by the transformation 


F (0) = P (l) (I„ + Q m (x )) F (1) . 


Thus, by Remark (7), Ga<i>(/ 1) has the form (C.5) with (C.6) and g = 1. 

By Remark 8, the matrix G A t n(/l = 0) is singular and the condition (C.7) does not occur, then one 
can find, by the same argument as above a permutation matrix and a nozero vector ur (2 '(x) in the 
left null space of G (1) (4 = 0). Let 


Q ( 2 ) (x) = [ < Ly > (x)] I <ij<n 


= 


s.t. 


D^lW 

\q 2) = 0 elsewhere 


for r + 1 < j <n — 1 


The matrix G A a,(A ) is then of the form (C.5) with (C.6) and g — 2. 

Consider the finite sequence of equivalent systems obtained by the transformation 


p(s-l) = />W(/ n + Qt s \x)) F {s) 


where 1 < 5 < g and 


= [q^(x)}\<uj<n 


s.t. 






-J+lJ 


(x) = -- 




„(*) 




= 0 elsewhere. 


for r+\<j<n-s+\ 


Then this process yields an equivalent matrix A(x, f A ) := A {,J) (x.f A < S )) with (C.6) for which either 
(C.7) occurs or g — n — r. But in the latter case one has, again by Remark 8, that det(Aj'(x)) = 0, 
and so (C.7) occurs. □ 


Example 12 (Continue Example 11). A simple calculation shows that det(Gx(d)) = 0 hence A 
is s-reducible. From (C.4), for A — 0, we have the singular matrix 


0 

0 

x + 1 

0 


1 

0 

0 

0 

x 2 

0 

0 

—2x 

. Let T — 

0 

1 

0 

0 

-X 

0 

0 

2 

0 

0 

0 

1 

0 

2 

0 

0 


0 

0 

1 

0 


then the transformation F — TG yields the equivalent system ff'ldG = A{x,f A )G where cr A = cr A 
and 


Mx,f A ) = 


& 

x- 2 

0 

—x 


-x 3 £a 

x£a 

2 

0 


0 

-2x& 
2 ft 


(1 + x)ft 

0 

0 

0 
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G a (A) has the form (C.5) with g = 1 and r — 2. In fact, 


G A (A) 


0 0 0 (1 + x) 

x 2 0 0 0 

0 2 A 0 

—x 0 2 A 


Lemma 8 . Given the system 


[AcrJ ? A x p dF = A(x, Za)F = Yj Ak(x ^ k A F - 

k =0 


Set r = rank(Ao(x )) and suppose that m E (A) > 1 and G a (A ) has the form (C.5) with conditions 
(C. 6 ) and (C.7) satisfied. Consider the shearing transformation S(f 4 ) = diag {fi A I r ,I n - r - e ,f A If) 
if q + 0 and S A ) = diag (f A I r , I n -r) otherwise. Then F — S (i; A )G yields the equivalent system 

CO 

[Act,] x p 8G - A(x, H A )C - 2 A k (x)f\ G, 

k =0 


for which rank(Ao(x )) < r and cr^ > 


(o- A if p > 1 
|cr A + otherwise. 


Proof. We partition A(x, g A ) as follows ( we drop (x, f A ) for clarity) 


A(x,£ a ) 


11 

A 12 

A 13 

21 

A 22 

A 23 

31 

A 32 

A 33 


where A 11 , A 22 , A 33 are of dimensions r,n — r — g , and g respectively. It is easy to verify then that 

A(x,£ 4 ) = S~ l AS-S~ l dS 

= S ~ l AS - % h A x p ~ l diag (crl r , O n - r - e , crl e ) 

A 11 -t A x p - x crl r ^A 12 A 13 

£4 A 21 A 22 I; a A 22 

A 31 ^A 32 A 33 - f h A x p ~ x cr A I e 


Hence, the new leading coefficient matrix is 


A 0 (x) = 


^0 

O 

M, 


U 1 O 

o o 

m 2 o 


and rank(Ao(x)) = rankiA^ Ui) < r. Moreover, due to the entries of the form ^ A x p 1 cr A in 
cr A should be adjusted as claimed if p < 1. 
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Algorithm 2 e-Rank Reduction of System [A :Ta ] for li > 1 
Input: h, p, cr A , A(x,^ a ) of [A (r J 

Output: R e GLJ'K) and an equivalent system [A^] which is e-irreducible. 
R <" In, 

h <— e-rank of A; 

f/(jc) <— Lemma 7 so that U~ 1 Aq{x)U has form (C.2); 

R <- /?[/; 

A <— U~ l AU - t^x p U~ l dU; Update cr A ; 
d = det(G^(A)); 
while af = 0 and /? > 0 do 
if h > 1 then 

T(x),g <— Proposition 4 ; 

A <- T X AT - % h A x p T- l dT\ Update cr A ; 

5(^a) Lemma 8; 

A <- 5“‘A5 Update cr A ; 

R<-RTS\ 

end if 

U(x) <— Lemma 7; Update cr A ; 

R <- RG; 

A <- U~ l AU - Z h A x p U- l dU\ 
d = det(G^(A)); 

/? <— e-rank of A; 

end while 
return (R, A,h). 
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Example 13 (Continue Example 12). Let S (= diag(£j,£j, 1, La) then G = S(£ A )IJ yields 
if h J)U = A(x,f)U where cr= = cr A and 


^4 

X" 

~ x £a 

0 

( 1 + x)e 

x h 

—2x 

0 

0 

2 h 


0 

— X 

0 

2 

0 


It is clear that the leading term Ao(x) has rank 1 < 2 = rank(Ao(x)). 

Proof. (Theorem 5) By Lemma 7, we can assume that Ao(x) is in the form (C.2). Then, 6 ) 4 (2) is 
constructed as in (C.3). Since det(GA(T)) = 0, it suffices to take the change of basis F = KG = 
TSG, where T and S are as in Propositions 4 and Lemma 8 respectively. □ 

Remark 9. The e-reducibility of[A aA \ implies that the rank of the leading coefficient matrix can 
be reduced (and consequently the e-Moser rank) without necessarily reducing the e-rank h of 
the system. If the e-reduction criterion holds for a sufficient number of equivalent systems then 
a repetitive application of such a transformation results in an equivalent system whose leading 
coefficient matrix has a zero rank, hence h can be reduced at least by one (e.g. Example 8). 


Appendix D. Proofs for Section 4 

We give the proof of Theorem 3 after establishing a series of useful lemmas. The following 
proofs are an adaptation to the parametrized setting of the proofs in [ 8 , Lemma 3, Lemma 4, 
Proposition 1, Theorem 1] respectively. Lor clarity within these intermediate proofs, we will 
express systems in the equivalent notation [A^J dF = A(x, s)F where A(x, e) e A4„(7C), rather 
than [Ao-J f h A x p dF = A(x,f A )F. Let 

det(T7 - A(x, e)) = A" + a„_i(x, e)/l„-i + • ■ • + ao(x, e). (D.l) 


such that a n = 1 and a,(x, e) = Yj°j= va i («■) ®u'W £; 6 7C for i e {0,We define the s- 
polygon N e (A) of [A^] as in Section 2.1.2, by taking P E (A) to be the union of P(i, val e (afx, e )) 
for i e {1, ..., n). We thus prove the following theorem, of which Theorem 3 is a straightforward 
corollary. 

Theorem 6. Consider the e-irreducible system [,4, ri ] dF = A(x, e)F where A(x, s) e Ai„(' K), 
h > 0, and (D.l). Ifh > n — rank(Ao(x )) then the s-formal exponential order is given by 


u> e (A) = max ( 

0 <i<n 


-val E (a,) 


)■ 


Additionally, the corresponding e-polynomial is given by the algebraic equation 


E S (X) = 


e 


Z 


a ik,val e (iXi t ) 


^(.k-io) 


where 0 < io < i\ <■■■< if = n denote the integers ifor which a> e (n - i) = —val e (ai ) (i.e. lie on 
the edge of slope o E of the e-polygon N E (A) o/tAo-J); and a Eva i t: ( ai )(x) — e^ vale (a,i ' afx, e)[ £ =o- 
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Lemma 9. Let A(x, e), W(x, e) be matrices in Ai n ('K) and Ai„CR) respectively. Put h — max((). —val e <A(x, e))) 
and let 

det(T7 - A(x, e)) - det(/l/ - A(x, s ) + W(x, s )) = a n - i(x, s)A "~ 1 + a„-2(x, s)A "~ 2 + ■ ■ ■ + ao(x, e), 
where a n -\(x, e), a n - 2 (x, s),..., ao(x, e) lie in *K. Then 


val K {a n -i) > (1 — i)h 1 < i < n. 


Proof. For any i e {1,..., n], it follows from Cramer’s rule that 

r i-i 

Otn-i = ^O/ Y\ a lJ 
1=0 s= 1 


where for 0 < l < y, 1 < .v < i - 1, w/ are entries in W(x,e) and a/ s are entries in W(x,s) or 
A(x, e). Consequently, 

val e {a n -i) > (i - 1 )val e (A(x, e)) > h. 


□ 


Lemma 10. Let A(x, e), B(x, e) be two matrices in A\„('K) such that val e (A(x, s)) < val e (B(x, s)). 
Consider the two systems [A^] dF — A(x,e)F and [TL-^] dG = B(x,e)G. Suppose that there 
exists T e GL n (fK ) such that [7L- S ] = T[A (Ta \. Put h = max( 0, -val e (A(x, e))) and let 

|det(/17 - A(x, e)) = A" + a„_i/l" _I + a„- 2 A"~ 2 + ■ ■ ■ + a?o 
|det(d/ - B(x, e)) = A" + /?„_i A n ~ l + /3 n - 2 A n ~ 2 + ■ ■ ■ + /? 0 . 


Then, we have: 

val E ((*n-i(x, s) -f3 n -i(x,s)) > (1 - i)h, 1 < i < n. 

Proof. By [31, Lemma 1] we can write T(x,e) = P(x,s) e y Q(x,s) where P,Qe GL n ('R) and 
e y — diag(e 71 ,... ,s y ") for some integers (y\ < yi < • • • < y„). Consider [A^] := PtA^] and 
[Ba- S ] := Q l [B aB \- Then we have 

B = e~ y As y , val e (A ) = val e (A), and valJB) = val e (B). 


It follows that, 

det(d/ - A) = det(/17 - P~ X AP + P l dP ) = det(A7 - A + (dP)P~ l ) 

■ det(d7 - B) = det(d7 - QBQ~ l + QdQ 1 ) = det(d7 -B + ( dQ l )Q ) 
det(d7 - B) - det(d7 - £~ y Ae y ) = det(/l7 - A). 

Since P, Q. P 1 , Q 1 , are units of Ai„CR), it follows that {dP)P~ 1 and (dQ l )Q inherit this prop¬ 
erty as well. The rest of the proof follows as a consequence of Lemma 9. □ 

Proposition 5. Consider the system [A, r J dF — A(x, e)F where A e M n ('K) and let: 

det(T7 - A(x, e)) - A" + ar„_i(x, £)A n ~ l + a„_ 2 (x, s)A"~ 2 H-+ a 0 (x, e). 
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Let [Cu- c ] dG = C(x, s)G where 


0 1 0 
0 0 1 


0 

0 


C(x,s) = Companion (c,(x, e))o<i<n-i = 


(D.2) 


0 0 0 

Co Cl C 2 


1 

c „—1 


be a companion system which is equivalent to [A^] over GL n (fK). Then we have, 
val e (oti — c,) > (1 — (n — i))h 0 < i < n — 1. 

Proof. Let [co-J with c n - 1 be the scalar equation representing [C, rc ] {cr c = crc )■ Consider k 
of [Co-J as dehned in Section 3.3. We define = (n - i)k, for i e {01} and [£>o- D ] = 
eP[ Co- c ]. Then D(x,s) = s~^C(x, s)eP where eP — diag(e®°,..., hP" 1 ). It follows that D(x,s) = 
s~ K Companion (fP'cf x, £))o<i<«-i • We have, val g (D(x, e)) ^ since val s (eP‘Ci ) = (n - z')/t + 
val(cj) > 0. By the equivalence between [A^] and [Co- c ] (resp. [Do- d ]) we have h > k (resp. 
val g (D(x, e )) > —k > -h). Let 

det(/17 - D(x, e)) = T" + d„-\(x, e)/l" _1 + ■ ■ ■ + e). 

Hence, by Lemma 10, we have 

(a, — dj) > (1 — (n - i))h, 0 < i < n. 

Moreover, by Lemma 9 

val E (di — t'i) > (1 — (n — z)) max (0, - val e (D(x,s))) > (1 — (n — if)K > (1 - in — i))h. 

It follows that 

val e (aj - Ci) > min (val e ( a - df), val E (dj - c,)) > (1 - (n - i))h. 

□ 

Proof. (Theorem 6) Let [ Ca - c ] dG = C(x, s)G be as in Proposition 5. Due to their equivalence, 
and [Co- c ] have the same e-exponential order and e-polynomial. Hence, we have 




n - i 


where 0 < z'o < z'i <•••< if = n denote the integer i for which a> E (n - i ) = -val E (ci). By 
Corollary 2 and Proposition 5 one has 


val E (cti - Cj) > (i — n + 1 )h — (i - n)u) E + (i - n)(h — a> E ) + h 

r 

> cj E (i - n) + (—n)(l-) + h > a> E (i — n) + r + h — n. 


n 


It follows from h + r — n > 0 that val E (a, - c,) > (i - n)u> E , 0 < i < n. 
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• If i e {i 0 ,... if) then aj e (i - n) — val e (c,) and val E (a, - c,) > val e (c,). Hence, val E (a,) — 
vciI E (C/) and (y l va / t ^ c E va i ^ c .y 

• Else val E (a, - c,) > oj E (i — n) and val E (c,) > to E (i — n ). Thus, 

val E (o',) > min ( val E (a, - c,), val E (c,)) > ai E (i - n). 


This completes the proof. 


□ 


We also illustrate Lemma 4 with the following example: 

Example 14 (Lemma 4). Consider the system [A (rA ] % 2 A dF = A(x, La)F with cr A - 0 and 


A(x,La) = 


0 0 L/C 0 U x 0 

L/C La 3 x%a 0 0 0 0 

0 0 0 x 0 0 4 

3 £tx 2 0 0 0 L/C xf A 0 

0 xf A 0 0 0 %a 2 0 

0 0 xU LaU 1 0 0 

0 La 0 0 0 0 xf A 


A(x,La) is e-irreducible, and we have n = 1, h = 2, and r := rank(Ao(x )) = 2. Thus, the 
condition h > n — r of Theorem 3 is not verified. Let us first consider a random ramification in La 
regardless of Lemma 4. For instance, let us consider La — La an> d a PPb' e-rank reduction. This 
yields system [Co- c ] P c dG = C(x,Lc)G with cr c — 0 and 


0 

0 

0 

xLc 

lc 

II 

0 

1 /4 L^xr 

x~l 2 c 

lc 

0 

1/4 Lc {P c ~ x 2 ) 

1/4 IJx 

3/4fc* : 

0 

0 

l| 

*l| 

0 

0 

II 

0 

4 

0 

0 

0 

0 

0 

*l| 

0 

0 

0 

0 

II 

31|* 2 

H 

0 

x~L 2 c 

0 

0 

0 

0 

X 

0 

0 

lc 

0 

II 

0 


for which, h = 5 and r — 2. Thus, the condition h > n — r is still not necessarily verified after 
a random ramification. To guarantee that we will arrive at a system verifying this condition, 
we make use of Lemma 4 and choose an integer d such that d > 7/(1 + 2/7) = 49/9. Let 
d — 7 > 49/9, perform La — P A [A<r A ]> an d then e-rank reduction. This yields the system 
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[fio-J £^'dG = B(x,f B )G with cr B — 0 and 


B(xJ b ) = 


& 

1/4 iv 1 

ps 

0 

i/ 4£ b (£^-* 2 ; 

1/4 ?£* 

3/4 ft 

0 

0 

0 

Ib* 

£8 

£8 

$B 

0 

0 

0 

£18 
? B 

IV 

0 

0 

£H 

4 

0 

0 

0 

0 

0 

0 

0 

IV 

0 

0 

0 

£H 


0 

£14 

iV 

0 

0 

0 

0 

0 

X 

0 

$B 

0 

£4 

0 

we have: h — 11, 

r = 2 

and so the condition h 

= 11 > n — 

r = 7-2 


so, Theorem 3 can now be applied. From det(/17 - we can compute 0 J vareps iioniM) = 

And so, we introduce the ramification f B = £ 3 and apply e-rank reduction which gives a s- 
system of order 42 and whose leading coefficient has 4 nonzero eigenvalues. In fact, we compute 
[M|^] £^<9W = M(x,^ m )W with cr M = 0 and 


M(x, | M ) 


14 


14 


28 


28 


|m * 

1/4 £ m x 2 

0 

0 

1/4 - l/4x 2 

1/4 |m * 

3/4x 3 





= 28 

s 28 


0 

0 

0 

X 



0 




s 14 



s 42 

0 

0 

0 

f,M X 

0 

0 


4 

0 

0 

0 

0 

0 

0 


= 28 




= 42 

z 14 , 

0 

* 

0 

0 

0 

|m 

3|m * 


s 14 




42 




14 


0 

0 


We remark that one can observe that this procedure was not necessary for such system since 
gcd( 14,42,28) = 14. Thus a simplification would leave us with a system of order 3. Thus, in the 
implementation we first try to ignore the condition h > n — r and apply Theorem 3 directly to the 
system [A^]. For this particular example we obtain a> e (A) = 3/2. Then by performing £4 = £^ 
and the e-rank reduction, we arrive at a e-irreducible system of order 3 whose leading coefficient 
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has 4 nonzero eigenvalues. In fact, we compute [A^] |^<9G = A(x,£)G where cr^ = 0 and 


AO,&) = 


0 

0 

0 

X 

P 

I 2 

0 

4 h* 2 

hx 

p 

h A 

0 

-l/4x 2 + 1/4 1\ 

1/4 x| 2 

3/4x 3 

0 

0 

P 

Ia* 

0 

0 

P 

b A 

0 

4 

0 

0 

0 

0 

0 


0 

0 

0 

0 

p 

^A 

3 h* 2 

3 

0 

Ia* 

0 

0 

0 

0 

X 

0 

0 

p 

b A 

0 

h 

0 


If Theorem 3 does not lead to the desired result, we then resort Lemma 4. In the case of unper¬ 
turbed system, it is also a matter of discussion whether this condition is necessary. M. Miyake 
recently claimed to have such an example (where this condition is necessary). 


Appendix E. Examples 

We treat with our algorithm examples from literature. For the computation of the full the 
e-exponential parts, we refer to our Maple package ParamInt [30], 

Example 15 (Continue Example 9). We resume the computation of the outer solution of system 
[Bo- b ] given by: 


[B rr J 


f B x 5 dW 1 = B(x,f B )W l = { 


0 1 


-1 

0 


1 -1 

0 0 

+ 

-1 

0 

% B + 

1 -1+x 4 


with cr B = -3. The leading matrix coefficient 
the characteristic polynomial of B gAj’ •' 


0 1 
0 0 


e B +o(f B )}w\ 

is nilpotent and s-irreducible. We compute 


l2 (o* 4 -o, (&- mV-\) 

A — --- A H— 


f B x 5 


Ax 10 


Then the s-formal exponential order is u> e = | and the e-polynomial is X 2 + ~ x =0. Thus, 
by setting e = e 2 , fs = f 2 B x 3 (fn = x~ 3 e) and applying e-rank reduction of Section 3.2 via 

diag(l,ffi), we get an e-irreducible system whose leading matrix coefficient is 

so, by the turning point algorithm of Section 3.1, we apply the transformation 


0 1 

-x -3 0 


. And 


1 0 

0 x- 3 ' 2 


which yields (note that we can also express the former using a change of independent variable 
x = t 2 ) 


[Bo-J P^x l9/2 dU = B(x,f B )U where cr B = -3 and 
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B(x,h) = 


0 1 


-x 3 / 2 0 

_ 

0 

0 


-1 0 

+ 

0 0 

\b + 

X 3 

0 

£~ + 


r 9/2 


0 


0 -*9/2 + *17/2 


| 3 « + 0(|1). 


The leading matrix has two distinct eigenvalues. Applying Splitting Lemma we get 
[Bo-,] f\x x9 ^ 2 dR — B(x, ftf)R where cr^--3and 


B(x,h) 


-i 0 
0 i 


_I r 3/2 

2 X 

0 


_ I j-3/2 
2 X 




Thus this system can be decoupled to any desired precision. The solution follows by straightfor¬ 
ward integration. Remark that, as expected, the system has the same s-formal exponential order 
and e-polynomials, as its equivalent scalar equation given in Example 2. The leading term of 

-2 i V* 


the e-exponential part Q is then given by: exp( E ^ 2 


0 

2 i yjx 


). Moreover, since crs. — -3, we 


set pi - 1/3. Further computations yield the full exponential part of outer solutions of the input 
system [Ao- a ] of Example 9: 

x 2 1 

—o- and 

2e- xe 

1 12 1 . 2 

- ^ g where x = -t . 

20 f 5 e 1 /“ 2tfe e 3 /- 

Moreover, since cr~ — —3, we set p\ — 1 /3. For the initial system of Example 9 given by: 


e z dF = 


0 1 0 

0 0 1 

s 0 x 


we have obtained so far outer formal solutions. To compute inner (and probably intermediate) 
solutions, we set r = xs~ PA = xe -1 / 3 which yields: 


s i d T F = 

which can be expressed equivalently as: 


0 1 0 
0 0 1 
E 0 T £3 


F, 






0 

1 

0 




\E<t e \ 

~f E d T F 

= E(T,h)F 

= 

0 

0 

1 

F, 

with cte = 





p 

\Se 

0 

T h. 




By Algorithm 1, 

one can 

compute the 

transformation F = 

TG 

where 




0 

0 

1 

(P 

$E 

0 

0 




T = 

0 

1 

0 

0 

is 

0 





1 

0 

0 

0 

0 

1 



?E- 


The resulting system is 

[E,rA 


iid T G = E(x,Ie)G = 

49 


r 0 1 
1 0 0 
0 1 0 


G, with erg = 0. 





































One can verify that Eq q has three distinct roots and so the system can be decoupled into three 
scalar equations. The leading term of the s-exponential part in the inner subdomain is then given 
by: 


f r dz+... 
£ 4/3 


exp( 


0 

0 


0 

f-l+i/ V3+...& 
2?73 

0 


0 

0 )• 
r-1-i/y/3+...dz 

2^ - 


And further computations show that the diagonal entries of the s-exponential part of a funda¬ 
mental matrix of formal solutions in the inner subdomain are given by: 


t + (l/ 6 )r 2 + (1/27)t 3 + 0(t 4 ) 

PV3 


and 


t RootOf(z 2 + z + 1) + (1 / 6 )t 2 - (l/27)(RootOf(z 2 + z + 1) + l)r 3 + 0(r 4 ) 

^ ' 


Since erg = 0, the reduction stops here. 


Example 16 (Iwano-Sibuya polygon). Consider the following scalar equation whose cr a = 0: 

KJ d 4 / - s-\2x + e 3 )d 2 / + s^x'f = 0, (E.l) 


One can verify that lo s = 2. The Iwano-Sibuya polygon is given by the following set of points: 
Paija) = (2,-1), Poo = (0, |), P 20 = (0, 5 ), ant/ P 23 = (§, 0 ) leading to only one slope given 
by cr = TTzen f/ze behavior can be investigated with the help of [p] 0 < pi = — 1 /cr = 

4/3. Hence, there is only one stretching transformation to consider: r = xe~ 4 / 3 . We now wz's/z 
to compute formal solutions and the sequence [p] wzY/z f/ze matricial representation using our 
algorithm: Let F — [f,s 2 df,s 4 d 2 f,s 6 d 3 f] T then (E.l) can be expressed as the following first 
order differential system (o~ A — 0) 


[A 0 


? A dF = A(x,Z a )F = 


1 

0 

0 


0 

1 

0 


0 + 2x 0 


F, where = s. 


(E.2) 


For the turning point treatment, we set x = t~ and then apply the transformation 


which yields the system 


T = 


0 

0 

,3 




fitdtG = A(t,f A )G = 


0 0 1 

0 t 2 0 

0 0 0 


0 0 


-3f 


4 + 2fit 

A 

2 

0 


A 

•3 *5/2 
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- 4 Q 

0 

0 


0 

0 

-23 

2 1 


0 

-2 1 
0 
0 


G. (E.3) 


















where cr^ — -3/2, i.e. £4 = t 3 ^ 2 e. 

Applying the Splitting lemma decouples (E.3) into two subsystems (truncated at order 4 in 
both t and s) 

[Sl sl ] gddtU 1 = S'fafoW 1 and [S* j2 ] f s 2 td t U 2 = S 2 (t,fa)U 2 
where cr s 1 = cr S 2 = -3/2, ij s i — 2 = t~ 2 ’i 2 s, and 


S\t,f s 0 = 


whose leading coefficient has two constant nonzero eigenvalues (hence the dimension of 
the system /order of the equivalent equation can be reduced). And, 

l/2f 2 sl (3f 4 sl t 2 +t 2 -4) r r 
1/2 t(2ft 2 t 2 + f 2 + 4) 72 ’ 


-1 /2f 2 sl (f 4 f 2 +t 2 + 6) 

-f 4 si t 2 + 2 

2f 2 f 5 ' 2 - f 4 f 2 - t 2 + 4 

-112^(3^?+t 2 + 8) 


S 2 (t,f S 2) 


(E.4) 


where 


fn = -i/s mfif 5 ' 2 - 9fy- - s^ 2 - 8) 


\r 2 = -l/16^ 2 (8^ 2? 5/ 2 _ io 5^ 2 -40^ 2 - 14f 2 - 16. 

Since et~ 2 ^ 2 — ex~ 2 i 4 , we have p\ = 4/3. 

If we wish to continue the reduction for the second subsystem, then, as usual, we first 
consider its leading coefficient: 


Slit) = 


0 t 

2t + 1/2 f 3 0 


Hence, by the turning point resolution (here it suffices to factorize t), the system [5 2 ] can 
be rewritten as: 

cA A a d _ 


tsy 

where crgi = - 2 , = r 2 e, and 

s 2 (t,f r -) = 

where 


17i =-1/23/ + 


9r£. 


^ 2 t A d t R = S\t,s)R 


l/2^ 2 (3l 4 4+l 2 -4) 7 ! 

t%_ + l/2t 2 + 2 72 


+%y 2 +1 


(72 = -l/16f 2 £? 2 (-105 t 4 f 4 2 + 8 f]/ - 40^ 2 f 2 - 14/ 2 - 16). 

The leading constant coefficient has two constant eigenvalues and hence can be decoupled. 

We consider again system [Ao- A ] which is given by (E.2), and we apply the stretching x — 
s 4 i 2 r which yields the system [iv E ]. One can verify that o~e = 0 and so with = e = 
= e 3 , we have: 


[£<rj 


t E d T H = 


51 


0 1 0 

0 0 1 

0 0 0 

-2r 3 ll 2 0 ~e E + 2r£ 


0 

0 

1 

E 0 


H. 


(E.5) 

















Then the s-rank reduction with 


yields 


[E r r e ] 


T = 


0 0 0 1 

0 0 || 0 

0 || 0 0 

II 0 0 0 


d T V = £(t,| £ )V = 


0 |l 5 + 2 r 0 

1 0 0 

0 1 0 

0 0 1 


.3 £12 


-2 T'c 

~ i 

0 

0 

0 


y 


°£ = 0,|e = S. 


Evidently, The e-exponential part oflE^] is zero. A formal fundamental matrix of solu¬ 
tions can be constructed following Subsection 2.1.1. 


Example 17 (Roo’s equation [33]). Consider the singularly-perturbed scalar differential equa¬ 
tion with cr a = 0 (f a — e): 

KJ ? a d 2 f = (x 5 + | fl x 2 + | 2 )/. (E.6) 


Iwano-Sibuya’s polygon o/fc/o-J consists of three segments connecting four points [33, p. 607], 
Poo = (0,5/2), P 0l = (1/2,1), P 0 2 = (1,0), one/ P„ e = (2, -1), and thus having the three slopes, 
—3, —2, and -1. Consequently, Iwano-Sibuya’s sequence [p] of positive rational numbers is: 


[p] 0 <pi = 1/3 <p 2 = 1/2 <p 3 = 1. 


We now wish to recompute the sequence [p] with the matrix representation o/[c/o-J using our 
algorithm. Let F = \f,f 2 df] T then [a (Fa ] can be rewritten as the following first order differential 
system (cr A =0) 


[VJ 


&dF=A(x,£ A )F = 


0 

x 5 +f A x 2 +f 2 A 



(E.7) 


• For the resolution of the turning point, let T — 
F — TG yields 


0 

x 5 ^ 2 


(resolution of turning point) then 




f 2 A x*dG=A(x,f A ) = 


1 

1+&+&X -\fi\x 5 ’ 2 


(E.8) 


with cr A = -3 andf A = x 3 e. The system [A (rA ] can be decoupled into two scalar equations 
using the Splitting lemma and we have p\ — 1 /3. 

We consider again the system [A aA ] and we perform the stretching r = xe -1 ^ and ramifi¬ 
cation /'/,■ = e = if — ||. This yields: 

[F ( r E ] ~f h: d T U = E(T,f E )U = 


52 


0 

t 5 i| + ||t 2 + h o 


U, with cr E = 0. 
















Let T i = _ g 3 (e-rank reduction) then U — T\G yields 


= E(t,%e)G - 


0 & 

ffi+T 2 0 


G, with cep = 0. 


The leading coefficient matrix is nilpotent. The computation of the e-exponential order 
suggests introducing an additional ramification in e of order 2. Since we already have a 
ramification of order 3, to avoid a proliferation of notations, we reset e — s 6 . This yields: 


qd T G = E(r,f~)G = 


0 4 

T 5 + T 2 + 0 

5 E 


with cr~ — 0, e — £e — s’* = e°. Let 73 

^ E 

s-irreducible system: 


= h 0 
o 1 


then G = TfW yields the following 


?aw= _5 , ; 2 


+ T- + fj 0 


W, with <jj — 0, and — s 6 — s. 


Let Tx — 


(resolution of turning point) then W = T^V yields 


ilrfi iyd T V = B(x,lj)V = r 3 + | + |2 _~y V 

with o~j — — 1, fj — As — T~ 1 g 1 ' ,6 . Clearly, system can be decoupled into two scalar 
equations. To find the restraining index, we rewrite r and s in terms of x and e: 

As = As 1 ' 6 = (x _1 £ 1/3 )£ 1/6 = x" 1 ^' 2 = (x~ 2 £) 1/2 , 

and so p 2 = 1/2. 

We consider again system [A^] and perform the stretching r = xe -1 / 2 , and the ramifica¬ 
tion s = s 2 . This yields: 


f c d T z = C(r,| c )z = 


AA + AA + £ 2 0 


with (Tc = 0, £(' - s = £ 2 = f 2 .. Lef T i 


n ?2 (e-rank reduction) then Z = T\G yields 

. ^cj 


: l/2 d T G = 


0 1 

tV /2 + t 2 + 1 0 


G, with ctq = 0 and - £c> 


which can be decoupled into two scalar equations by Splitting lemma. Since <Jq = 0, the 
restraining index is infinity. 

And so, by our techniques, we have obtained the two slopes pi = 1/3 and P 2 — 1/2. 

53 



Example 18 (([16], Example 1, Section 9.5, p. 446)). Consider the following singularly- 
perturbed scalar differential equation 


ed 3 f -df + xf = 0 

which can be rewritten as the following differential first order system 


[A<rJ 


€a8F = A{x,S a )F = 


0 

0 

-x 


if a 0 

0 £t 
1 0 


where F — [f,df,d 2 f] T , cta = 0, and £4 = e. Since Aq(x) is nilpotent, Splitting lemma 
cannot be applied but no treatment of turning points is required. The transformation F — 


G yields the following e-irreducible equivalent system: 


0 

1 

0 


1 

X 

0 

1 

0 

0 


0 

1 

0 

0 

0 

1 


0 

0 

1 




f A dG=A(x,f A )G = 


-£?* 

£4 

1 


£4 - £? & 

f A x 0 

0 0 


G, 


with cr A — 0 and 8 A (A) = —A - x. Since the leading matrix coefficient is still nilpotent, we have 
to compute the e-exponential oder. The characteristic polynomial of is given by: 


A 3 + ^ ~ ' A + 


& 


X 

Ia 


Consequently, the e-exponential oder is given by: a> B (A) — max {^, 5 } = Let e = e 2 and 
f — f 2 then we have: 


-f l x -x 2 | 2 -^ if 


£.2 


[A^] fdG 

= A(x,fi)G = 


| 0 0 

0 

The transformation G = 

0 | 0 

1 


0 0 1 

0 


system: 


I 2 

| 2 x 

1 

0 

1/2 

1/2 

0 

0 

- 1/2 

- 1/2 . 


0 

0 


G, with cr * = 0. 


U yields the following e-irreducible 




fax 

hi 2 

h/2 

[B ( r„] 

f B dU = B(x, £h)U = 

{-x 2 -\)h 

- 1/2£bx — 1 

-1/2 l B x 



(— x 2 - i)£b 

— 1/2£b x 

— 1/2£b x + 1 
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with cr B — 0, fo — s. The leading matrix coefficient Z?o,o has three distinct roots and consequently, 
the system can be decoupled into three first order scalar differential equations by applying the 
Splitting lemma. A fundamental matrix of formal solutions is then given by: 


i> 0 , s 1/2 ) exp( 


0 0 0 

0^0 

0 0 -fc 


where <I>(x, e 1 ^ 2 ) is the product of all the transformations performed including that of the Splitting 
lemma. Since cr B = 0, the origin is not a turning point for this system and we do not need 
stretchings. 

Example 19 ([? ]). We consider the scalar differential equation 

M s 3 d 2 f - (x 3 + e)f = 0 , 

with cr a = 0. Setting F = [/, df] T , we get the following differential first order system 




? a 8F = A(x,f A )F = 


0 

x 3 + fo 0 


with cr A — 0 and fo — s. 


Since Aqq is nilpotent and Ao(x) is not, we first start with the treatment of the turning 
point at x — 0. It suffices however to factorize x 3 from Aq(x) which results in cr A = -3 


and fo — x 3 s. We then apply the transformation F — 
system: 




£x b dG = A(x,fo)G = 


fo 0 

0 1 


G to get the s-irreducible 


3x 3 fo fox 0 
1+fo 0 


with o~a = _ 3 an d ~ x is - Since Ao(x) is still nilpotent, we proceed to computing the 
s-exponential order from the characteristic polynomial of which is given by: 

x fox 6 


Hence, cl> £ (A) = j and the s-polynomial is given by: X l —x i = 0. Let s = ef and ^ = x 5 s, 
then with fo = fox 3 we have: 




x 12 fo<9G = A(x,fo)G = 


1 + x 3 fo 0 

^A 


with o~a — —3. We then perform G — 


fo 0 

0 1 


U to get the e-irreducible system: 


[Bcr B ] 


x n f B 8U = B(xJb)U = 


3x"| 2 (fo + 1) 
l+x 3 f 2 


0 


u. 


with cr B = -3 and fo = x 3 e. 
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Now that Bq(x ) has two distinct eigenvalues and Bop is nilpotent, we perform U = 
to treat the turning point. Then we have: 


x 9 ' 2 

0 


0 

1 


w 


[Bcr J 


x 15/2 AdW = B(x, £~h)W = 


t-x% 0 


1 


W, 


with crft — -3 and gg = x 3 e. 

Since Boo has two distinct eigenvalues, ± 1, one can proceed to the Splitting lemma which 
decouples the system thoroughly. From the the eigenvalues of Boo or the E-polynomial, one 
can read the e-exponential part of an outer solution of this system: 


exp( 


fz ,/2 dz 

e 3/2 

0 


fz ,/2 dz 

s 3/2 . 


and deduce that pi = 1 /3. 

• We consider again system [A o- A ] and perform the stretching r = xs~ [ ^. After introducing 
the ramification e — e 3 , we have: 


[Soj] 


^d T U = S(T,fs) = 



it 

0 


U 


with cr$ = 0, fs = s, fs = Is- The transformation U - 


p 

0 


0 

1 


V yields 


[^ s ] 


l 2 d T V = S(T,is)V = 


0 

T 3 + 1 


1 

0 


y 


with erg = 0, £g = s. Due to the nature of the eigenvalues of Sop, the system [S^-,] can be 
decoupled into two scalar equations, erg — 0 suggests that there are no more stretchings 
to perform. 


However, if we consider again system [A^] and experiment with the stretching r = xs *, 
we get: 

0 #| 

?£ + i o 


[Soj] 


£sd T U = S(j,£s)U = 


U 


with cr s = 0 and fs — £ - By the transformation U = 
get 

[So-,] d T G = S(T,fg)V = 


E 0 
0 1 


V and some simplification we 


0 


T 3 ^ + 1 


y 


with erg = 0 and fg = e. Due to the nature of the eigenvalues of Sop, the system [5o- s ] can 
be decoupled into two scalar equations. Its E-exponential part is zero. 


So far, this matricial approach can determine the stretching r = xe -1 ' 3 but not the stretching r = 
xe~ 1 . However, both stretchings can be obtained from Iwano-Sibuya’s polygon by the treatment 
of this system as a scalar equation [1 , Last section]. 

56 



Example 20 ([16], Example 1, Section 9.5, p. 453). We consider the scalar equation 


e'xd 2 f + x 2 df - (x 3 + s)f = 0, 

which can be rewritten as the following first order differential system 


[ArJ 


xf 3 A dF=A(x,f A )F = 


0 




X'’ + f A - X - 


where cr A = 0, f A = s, and F = [/, <9/] r . We first start with the treatment of the turning point at 
x = 0 since Aq.o is nilpotent while Aq(x) is not. It suffices here to factorize x 2 which yields 


0 £ 3 x 5 

x + & -1 


[Ao-J x f/IF = A(x,f A )F = 

where cr A — -2 and f A — x~ 2 e. Now that Aq.o has two distinct eigenvalues, we can proceed to 


the Splitting lemma. The transformation F 


0 1 
1 x 


G block-diagonalizes Aq(x) which yields 


lArJ 


x 5 f]dG = A(x,if~)G = 


A 

_ v -2 


-fix 6 -! * 


3 v 5 


3* 


A ,4 

f]x 6 


G, 


where cr^ = -2 ant/ = x z s and from which we can deduce that the e-exponential part is 
given by: 

ffr + 0(ln(x)s~ 2 ) 0 

0 0(ln(x)£~ 2 ) 


exp( 


)■ 


Since the restraining index is nonzero, we can deduce that p\ = 1/2 and apply the stretching 
t = xe ^^ 2 in the input system [A^] and repeat the formal reduction procedure. 

Example 21. In this example, we use our algorithm to explain some of the computations in the 
introductory example in the light of the eigenvalues of the matrix of the system. Thus, we consider 
again the system [A^] given by (1); 


[ArJ 


s dF = A(x, s)F 


; 

X - £ 


F. 


We remark that the eigenvalues of A(x, e) are given by /l 12 = ±(x 3 - e) 1 ^ 2 . We thus have the 
two possible expansions which correspond to the two subdomains above: 

• 4 12 = ±(x 3 - e) 1/2 = +/e 1/2 (l - ^£ _1 x 3 - 4e _2 x 6 + ...), valid for |£t'x 3 | < 1 . 


A \2 = ±(x 3 - e) 1/2 = ±x 3/2 (l - \ex 3 - ie 2 x 6 + ...) , valid for \ ex 3 [ < 1 . 
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